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Abstract. Price's Law states that linear perturbations of a Schwarzschild black hole fall off as 
t~ for t — y oo provided the initial data decay sufficiently fast at spatial infinity. Moreover, if 
the perturbations are initially static (i.e., their time derivative is zero), then the decay is predicted 
to be t~ 2e ~ i . We give a proof of t~ 2e ~ 2 decay for general data in the form of weighted L 1 to 
L°° bounds for solutions of the Regge-Wheeler equation. For initially static perturbations we 
obtain t~ 2e ~ 3 . The proof is based on an integral representation of the solution which follows from 
self-adjoint spectral theory. We apply two different perturbative arguments in order to construct 
the corresponding spectral measure and the decay bounds are obtained by appropriate oscillatory 
I '— ^ , integral estimates. 
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1. Introduction and main result 

O . 

In General Relativity, the dynamics of spacetime is governed by Einstein's equation which, in 
the absence of matter, takes the form 

R^u(g) = o 

where R^ v {g) is the Ricci tensor of the Lorentz metric g. Exact solutions (i.e., solutions which 
are known in closed form) include the free flat Minkowski spacetime as well as the Schwarzschild 
(NJ ' metric and, more generally, the Kerr solution. The Schwarzschild solution is spherically symmetric 

and corresponds to a nonrotating black hole whereas rotating black holes are described by the 
axially symmetric Kerr spacetime. A fundamental mathematical problem in General Relativity is 
the understanding of the stability of these solutions. The stability of the flat Minkowski spacetime 
under small perturbations was shown in the seminal work of Christodoulou and Klainerman [15] in 
the late 1980's. A simpler proof was later developed by Lindblad and Rodnianski [37]. However, 
we are very far from understanding the dynamics near a black hole. Yet, latest experimental 
setups are crucially dependent on such an analysis, in order to observe gravitational waves (see 
for example [21], [23], [23], [22] and cited ref.). Most efforts are now focused on understanding 
the linear dynamics and stability of such solutions, see e.g. [32j, [7] and cited ref., as well as [33]. 
The mathematical aspects of the problem will be discussed below in more detail. We also refer the 
reader to the survey [20J which gives an excellent overview of recent developments in the field from 
the mathematical perspective. 

1.1. Wave evolution on the Schwarzschild manifold. As a first approximation to the lin- 
ear stability problem of a nonrotating black hole one may consider the wave equation on a fixed 
Schwarzschild background. One is then typically interested in decay estimates for the evolution. To 
simplify things even more, one restricts the analysis to the exterior region of the black hole which, 
however, is physically reasonable: such a model describes a black hole subject to a small external 
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perturbation by a scalar field — a situation which, with a more realistic matter model, is certainly 
relevant in an astrophysical context. In order to formulate the problem we choose coordinates such 
that the exterior region of the black hole can be written as (t,r, (9,<ft)) £ R x (2M, oo) x S 2 with 
the metric 

g = -F(r)dt 2 + F(r)~ l dr 2 + r 2 (d9 2 + sin 2 9d(j) 2 ) 

where F(r) = 1 — and, as usual, M > denotes the mass. We now introduce the well-known 
Regge-Wheeler tortoise coordinate r* which (up to an additive constant) is defined by the relation 

F = — . 

dr* 

In this new coordinate system, the outer region is described by (t, r*, (9, 0)) G R x R x S 12 , 

(1) 5 = -F(r)dt 2 + F(r)dr 2 + r 2 (d9 2 + sin 2 (9# 2 ) 

with F as above and r is now interpreted as a function of r* . Explicitly, r* is computed as 

r* = r + 2Mlog(^-l) . 

Generally, the Laplace-Beltrami operator on a manifold with metric g is given by 

and thus, for the metric g in (pQ), we obtain 




By setting i/)(t,r*,9,(f)) = r(r*)i/>(i, r*, 0, 0) and writing x = r*, the wave equation = is 
equivalent to 

(2) - d 2 ^ + a 2 ^ - -—V + -A 52 V = 0. 

r ar r z 

The mathematically rigorous analysis of this equation has been initiated by Wald [58], however, 
the first complete proof of uniform boundedness of solutions is due to Kay and Wald [34] . Recently, 
Dafermos and Rodnianski have found a more robust method to prove boundedness of solutions 
based on vector field multipliers that capture the so-called red-shift effect [15] ; see also [20] for 
a survey and generalizations of these results. The goal of our present work is to prove L 1 to L°° 
decay estimates for Eq. ([2]). Different types of decay estimates have been proved before. Local 
decay estimates, based on multipliers generalizing the Morawetz estimates, were initiated in [SJ, [9] 
and [ID]. Later, a similar approach was used in [IT], [IT], [18], [19] and [38] to prove both local decay 
estimates and pointwise decay in time based on conformal type identities. In [12], [41] , |55j and [39 
it is shown how to apply such estimates to obtain Strichartz type decay estimates. We also mention 
the recent work [3]. After submission of the present paper, Tataru announced a proof of the sharp 
pointwise t -3 decay for general data without symmetry assumptions, see [54] , Moreover, his result 
also applies to the more complicated case of rotating Kerr black holes. In fact, in the follow-up 
paper [27] we also obtain pointwise t~ 3 decay on Schwarzschild for general data. We will discuss this 
below in more detail. Our results differ from the above in certain respects: the methods we use are 
based on constructing the Green's function and deriving the needed estimates on it. Previous works 
in this direction include mainly the series of papers [29] . |30j . [28] where the first pointwise decay 
result for Kerr black holes has been proved, see also [35] and [36J for Schwarzschild. In our approach, 
we freeze the angular momentum i or, in other words, we project onto a spherical harmonic. More 
precisely, let Y^ m be a spherical harmonic (that is, an eigenfunction of the Laplacian on S 2 with 
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eigenvalue —£(£ + 1)) and insert the ansatz ip(t, x, 9, (ft) = ifti, m (t, x)Yi tTn (9, eft) in Eq. ([2]). This yields 
the Regge-Wheeler equation 

d?iftt,m ~ d%.ipi,m + Vi ia (x)ifti, m = 

with a = 1 where 

2M\ (£(£ + 1) 2Ma\ 
r(x) J \ r 2 (x) r 3 (x) / 

is known as the Regge-Wheeler potential In the present work, we obtain decay estimates for 
solutions of this equation. However, before we explain our results in more detail, we further motivate 
the study of the Regge-Wheeler equation by considering more general black hole perturbations. 

1.2. Other types of black hole perturbations. The wave equation on the Schwarzschild man- 
ifold describes the time evolution of linearized scalar field perturbations of a black hole. Of course, 
not all physically relevant situations are covered by this simple model since it ignores the underlying 
tensorial structure altogether. Eventually, one is interested in perturbing fields of higher spin, in 
particular gravitational perturbations. However, as a remarkable fact, the Regge-Wheeler equation 
is also relevant in this context. This follows from a reduction procedure that goes back to Regge 
and Wheeler [49J as well as Zerilli [59] , see also [57] and [2] . We will briefly sketch how this comes 
about. In order to study gravitational perturbations, one considers a perturbed Schwarzschild 
metric g of the form 

g = _ e 2{u+5u) dt 2 + e 2(^+^) (# _ 6udt _ dg2dr _ 5qzd9) 2 + e 2^ 2 +S^) dr 2 + ^3+^3)^2 

where the various coefficients are allowed to depend on t, r, 6 and e 2u = e~ 2 ^ 2 = 1 — e^ 3 = r, 
= r sin 9 (we follow the notation of [14]). It can be shown (see [13]) that this ansatz is sufficiently 
general. One then requires the metric g to satisfy the linearized Einstein vacuum equations, i.e., 
one linearizes R^ u {g) = with respect to the perturbations 8v, 5tp, etc. It turns out that one 
has to distinguish between so-called axial (5co,5q2,Sq3) and polar (5v,5ip,5iJ,\,5n2) perturbations, 
depending on the behavior of the metric under the reflection (ft i— > —(ft. After a lengthy calculation 
and separation of the ^-dependence one arrives at 

£(£ + 1) 6M\ _ fl 

r 2 (x) r 3 (x) / 

where ipi is an auxiliary function which completely determines the axial perturbations, see |14| 
for details. Thus, ifti satisfies the Regge-Wheeler equation with a = —3. In the case of polar 
perturbations, Zerilli [59] has derived an analogous equation with a more complicated effective 
potential. However, Chandrasekhar [13] (see also [13]) has found a transformation involving differ- 
ential operations that relates this equation to the one for axial perturbations. As a consequence, 
the Regge-Wheeler equation provides a fairly complete description of gravitational perturbations 
with a fixed angular momentum parameter £. Moreover, we mention the fact that the Regge- 
Wheeler equation with parameter a = appears in the study of electromagnetic perturbations of 
Schwarzschild black holes, i.e., if one considers the Einstein-Maxwell system and linearizes around 
the Reissner-Nordstrom solution with zero charge. We do not comment on this further but simply 
refer the reader to the literature, see [14] and references therein. As a consequence, the study of the 
Regge-Wheeler equation can provide valuable information on the stability of Schwarzschild black 
holes under various types of perturbations and it is truly remarkable that such a unified approach 
is available. 
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1.3. Decay estimates for the Regge Wheeler equation. The most salient feature of the 
Regge-Wheeler potential is that it decays exponentially as x — > — oo which corresponds to ap- 
proaching the black hole, whereas for x — > oo, it falls off as x~ 2 . Strictly speaking, this is only true 
for I > 0. The case I = is exceptional and we consider it separately in the companion paper |26| 
where we obtain the sharp i~ 3 decay as predicted by Price's Law. Consequently, in this paper, we 
focus on £ > which, unless otherwise stated, will be assumed throughout. For £ > the potential 
has inverse square decay and it is well-known [25J that this fall-off behavior is in some sense critical 
for the scattering theory. In order to explain this we define the Schrodinger operator CT by 

He, a f = -f" + Vt, f 

and recall that the Jost solutions /±(x, A) are defined by % <r/±(-,A) = A 2 f±(-, A) and f±(x,X) ~ 
e ±i\x ag x _j. -j- 00 _ xhe property Ve^ a G L X (IR) is sufficient to guarantee the existence of these 
solutions, see [25], but the inverse square decay of Vg^ a is critical in the sense that at this power 
the Jost solutions typically are no longer continuous as A — > 0. Nevertheless, following [53], it is 
possible to perform a detailed spectral and scattering analysis of the Schrodinger operator / H^ (J . 
However, we emphasize that the present work differs considerably from [53] due to the asymmetric 
decay properties of the potential V^ CT . Of particular importance is the asymptotic behavior of the 
resolvent ((A -MO) — %£ )(T ) _1 (and thus, of the Jost solutions and their Wronskian) as A — > 0. This is 
a common feature in dispersive estimates, see [51]. In particular, we are faced with the possibility 
of a zero energy resonance. However, it was already observed earlier [35], [16], [33], [57] that in 
the physically relevant cases such a zero energy resonance does not occur (see also Lemma 16.31 
below). Our approach is detailed enough to show rigorously, for the first time, the decay estimates 
depending on the angular momentum of the initial data. In his seminal work [33], [35], see also 
[31] . [35] , Price heuristically derived the decay rate in time at a fixed point in space, and concluded 
that, depending on initial conditions, the decay rate is either t~ 2e ~ 3 or t~ 2e ~ 2 where £ is the angular 
momentum. This result is now referred to as Price's Law. There has been some confusion in the 
literature concerning the precise prediction of Price's Law. This has been clarified in Price and 
Burko |46j . If the initial data decay sufficiently fast at spatial infinity then the pointwise decay 
in time is predicted to be t~ 2i ~ 3 . In the present paper we give the first proof of an ^-dependent 
decay rate. More precisely, we obtain a t~ 2i ~ 2 estimate which is one power off the sharp version of 
Price's Law. However, we emphasize that our method yields estimates in terms of the initial data 
and not just a pointwise decay law as is common in the physics literature. To be more precise, we 
show that 
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for all t > where w a (x) := (x) a is a polynomial weight and, as usual, (x) := (1 + \x\ 
Here one has to require 1 < a < 2£ + 3 and one needs to exclude Q (£,a) G {(0, 0), (0, -3), (1, -3)} 
which are exactly those cases where zero energy resonances do occur — however, they are physically 
irrelevant due to a gauge freedom, cf. [14]. Observe that for a = 2£ + 3 we obtain precisely the 
aforementioned bound. It is also obvious from our approach that the decay of initially static 
perturbations is better by one power of t as is reflected by our cosine estimate. This is a general 
effect which is also present in Price's prediction, cf. [46] . 



In fact, the case (£, a) = (0, 0) is trivial since for these parameter values the Regge-Wheeler equation reduces to 
the free wave equation on the line. 
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The proof of ([3|) and @ is based on representing the solution as an oscillatory integral in the 
energy variable A, schematically one may write 



rp(t,x) = J U(t,X)lm [G e ,a(x, x' , A)] f(x')dx'd\ 

where U (t, A) is a combination of cos(tA) and sin(tA) terms and G^ a {x, x' , A) is the kernel (Green's 
function) of the resolvent of the operator / H^ IJ . Gi j(7 (x, x' , A) is constructed in terms of the Jost 
solutions and we obtain these functions in various domains of the (x,X) plane by perturbative 
arguments: for |xA| small we perturb in A around A = 0, whereas for |xA| large we perturb off of 
Hankel functions. This is done in such a way that there remains a small window where the two 
different perturbative solutions can be glued together. One of the main technical difficulties of the 
proof lies with the fact that we need good estimates for arbitrary derivatives of the perturbative 
solutions. This is necessary in order to control the oscillatory integrals. The most important 
contributions come from A ~ and we therefore need to derive the exact asymptotics of the 
Green's function and its derivatives in the limit A — > 0. For instance, we prove that 

Im [G v (0, 0, A)] = XP e (X 2 ) + 0(X 2i+1 ) 

as A — > 0+ where Pi is a polynomial of degree t—\ (we set Pq = 0) and the O-term satisfies 
0^ k \X 2i+1 ) = 0{X 2t+1 ~ k ) for all k G No- Our approach therefore yields further information on the 
Green's function and the fundamental solution of the wave equation on the Schwarzschild manifold. 



1.4. Interpretation of the result and further comments. For the relevant parameter values, 
i.e., (£, a) {(0, 0), (0, — 3), (1, — 3)}, the Regge- Wheeler potential Vt^ is positive, decays as de- 
scribed above, and has a unique nondegenerate maximum at r = r$ which is known as the photon 
sphere. As a helpful analogy, consider a Newtonian particle placed at r = tq with vanishing kinetic 
energy but potential energy V^ jO -| r=r0 . It will remain at rest, but every slight perturbation will 
make it lose potential energy and gain kinetic energy; the larger £ is, the faster this will occur 
which reflects itself in the ^-dependent decay rates. In the context of the Schwarzschild geometry, 
the flow of null geodesies near r = r$ is unstable and the dispersion provides a mechanism that 
spreads out the wave away from the photon sphere. Moreover, the higher the angular momentum 
I, the faster the geodesies will pull away leading to the accelerated decay provided by Price's Law. 
Note carefully, however, that this is counteracted by what can be viewed as a stabilizing effect of 
large I. Technically speaking, this reflects itself in the constant C£ tCC : the larger this constant is, 
the longer one has to wait before the decay estimates become effective. It is important to note 
that our approach is essentially blind to the local geometry, that is, the fine structure of the po- 
tential is irrelevant — only positivity, the decay properties and the nonexistence of a zero energy 
resonance are used. This is in contrast to the methods based on Morawetz type estimates. In 
particular, the phenomenon of trapping does not play a role at this level — it simply enlarges the 
constants. However, eventually one is interested in the overall decay which is obtained by summing 
the individual contributions over all £ and at this stage, of course, trapping becomes relevant since 
the ^-dependence of the constants is crucial for the summation. As a matter of fact, our proof 
produces a constant which grows super-exponentially in I. Consequently, in order to be able to 
sum the estimates, a different approach is necessary for large i. This issue is addressed in our paper 
|27j where a detailed semiclassical asymptotic analysis is performed. The role of the semiclassical 
parameter h is played by £ _1 (simply divide He a by I 2 ). In particular, such an analysis requires a 
careful study of the spectral measure near the maximum of Vi a and it is exactly at this point where 
the instability of null geodesies at the photon sphere becomes crucial. As a consequence, in [27], we 
show that the estimates for individual €s can indeed be summed and thereby, we obtain the sharp 
i~ 3 decay bound for general data with a loss of a finite number of angular derivatives. We also 
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remark that there are various formal approaches in the physics literature to find the ^-dependence 
of the constants, see in particular [3]. 

Let us finally mention that decay estimates like (J3j) and @ play an important role in the current 
theoretical and numerical analysis of black holes. For instance, they serve as a way to verify various 
numerical schemes for solving the Einstein equations in the presence of black holes, see e.g. [47J, 
0, @H], [6], 0, [3, [12], @], [21], [23], [21], [22] and cited ref. For other recent theoretical 
implications of the angular behavior see for example [I], [6] and [40] . 

1.5. Notations and conventions. For a given smooth function / we denote by 0(f(x)) a generic 
real-valued function that satisfies \0(f(x))\ < \f(x)\ in a specified range of x which follows from 
the context. We write Oc(f(x)) if the function attains complex values. The symbol f(x) ~ g{x) 
for x — > a, where g is smooth, means lim x ^ a = 1. Furthermore, the letter C (possibly with 
indices) denotes a generic positive constant. We say that 0(x 7 ), 76I, behaves like a symbol, is of 
symbol type, or has symbol character, if the k-th derivative satisfies 0^ k '{x^) = 0(x 7_fc ). As usual, 
we use the abbreviation (x) := \f\ + \x\ 2 and the symbol A < B means that there exists a C > 
such that A < CB. We also note that all of the functions we are going to construct will depend 
on the parameters £ and a, however, in order to increase readability of the equations, we will omit 
this dependence in the notation most of the time. The same comment applies to all implicit and 
explicit constants. Finally, as already mentioned, we assume I > unless otherwise stated. 



2. Solutions of the Regge- Wheeler equation 

2.1. Asymptotics of the potential. As explained in the introduction, linear perturbations of 
the Schwarzschild spacetime are described by the Regge-Wheeler equation 

iptt ~ tpxx + V£ )(T (x)i/j = 

where V^ CT is the Regge-Wheeler potential and x is the tortoise coordinate which is related to the 
standard r-coordinate by 

(5) x = r + 2Mlo g y^-l) . 

Vi t(J is given by 

v ,, = (i_ + + 2M ° \ 

,a \ r(x) J \ r 2 (x) r 3 (x) / 

where r{x) is implicitly defined by Eq. ([5]). The valid range of the parameters is I G No and 
a = —3, 0, 1. We start by obtaining the asymptotics of the potential V^ a . 

Lemma 2.1. The function x 1— > r(x) has the asymptotic behavior r{x) = x(l + 0(x~ 1+e )) for 
x^oo and r(x) = 2M + 0{e x / {2M ^) for x — > —00 where e G (0, 1) is arbitrary and the O-term in 
the expression for x — > 00 behaves like a symbol. 

Proof. The function r{x) is implicitly defined by x = r{x) + 2M log (^jj- — 1^ and thus, we have 
r(x) —7- 00 as x —7- 00. This implies x ~ r(x) and hence, r{x) ~ x as x — > 00. We infer that 
x - r(x) = 2Mlog - l) ~ 2Mlog ^ ~ 2Mlog gfy and this shows x - r(x) = 0(x £ ). For 

the symbol behavior note that ^p(r) = (l — 1 which implies that r'(x) = 1 — ^pj- The claim 
now follows by induction. 

For the case i^-oowe have e x ^ 2M ^ = e r ^/ (2M \^--l) ~ e(^-l) which shows r(x)-2M ~ 
2Me x /( 2M ^>~ 1 and this implies the claim. □ 
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Corollary 2.1. The Regge-Wheeler potential Vf >cr has the asymptotic behavior 

^(x) = ^^(l + 0(x- 1+£ )) 



ar 

/or x — t- oo and V£ jCr (x) = 0(e x '( 2M ') for x — > — oo where e G (0, 1) is arbitrary and the O-term in 
the expression for x — > oo behaves like a symbol. 

Proof. Just insert the asymptotic expansions from Lemma 12.11 in the expression for V^ i(T . For the 
symbol behavior apply the Leibniz rule and Lemma IA.11 □ 

2.2. Hilbert space formulation. We define the Schrodinger operator Hi j(T on L 2 (M) with domain 
V(H e>a ) := H 2 (R) by 

Hi, a f := -/" + V itV f. 

From the decay properties of Vi^ it follows that is self-adjoint (see e.g. [56]). Furthermore, 
integration by parts shows 

{U^f\f) L 2 m > (V^f\f) Lm > 

since Vg iCr > for all £ G No if o - = 0, 1. For gravitational perturbations (o~ = —3) we have to assume 
£ > 2 to obtain Vt^ a > which we shall do from now on. We conclude that the spectrum of %t jCr 
is purely absolutely continuous and we have a(7ii tCT ) = a ac (T-li i(T ) = [0, oo) (see [56]) provided that 
(£,cr) ^ {(0, —3), (1, —3)}. An operator formulation of the Regge-Wheeler equation is given by 

where VP : M — > L 2 (R). Applying the functional calculus for self-adjoint operators, the solution \& 
with initial data ^(0) = / and ^jr(O) = g is given by 



nt) = co S (t v m^)f + Mt ^ ) 9. 

V Ttl,a 



Thus, in order to obtain decay estimates for the solution, we have to understand the operators 

S'm(ty/H£ }a ) 



cos(ty/TLi^) and 



v 7 ^ 



2.3. The spectral measure. Recall that the spectral theorem for self-adjoint operators asserts 
the existence of finite complex-valued Borel measures fj, u ^ v such that, for u,»£ T>(7ii a ), we have 



{Hi,aU\v) L 2(u) = \ Xdfi UiV (X). 
J 

The solution operator cos(ty / ^ jCr ) is then given by 

i(ty/rli ta )u v) = I cos(tv / A)o!/i U)1 ,(A) 



cos( 

for u,v G L 2 (M) and analogous for the sine evolution. The point is that the spectral measure can 
be calculated in terms of the resolvent R% lcr (z) = (z — T-L^a) OI Indeed, for u G L 2 (M) set 

f°° 1 

F u {z) := -(R Hea (z)u\u) L 2 {R) = / dn u {\) 

Jo <\ — z 

where [i u := fj, u>u and Imz > 0. F u is the Borel transform of the measure \i u and, since the measure 
\x u is purely absolutely continuous, we have 

dfJ, u (X) = — lim Im(F u (X + ie))dX, 

7T e— >0+ 
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see [56] for the underlying theory of this. The measure fj, U)V can be reconstructed from \x u by the 
polarization identity, i.e., fi UjV = ^(fi u+v — [i u -v + i^u-iv — i^u+iv)- Furthermore, the resolvent is 
given by 

R Ul a {z)u{x) = / Gi, a {x,x' ,y/z)u(x')dx' 



where Gi a is the Green's function (we always choose the branch of the square root with Im-^/z > 
if Imz > 0) and thus, we have 



d/i u (X) 



- — lim 

7T £->0+ 



Im 



G^ )0 -(x, x', VX + ie)u(x')u(x) 



dx'dxdX. 



It is known (and, for the convenience of the reader, will be shown below) that the limit 

Gi a (x, x', vA) := lim Gi CT (x, x' , y/X + ie) 

exists and satisfies sup^ x ; g]R \Ge >a (x, x' , A)| < 1 for all A > Ao where Ao > is arbitrary. Thus, if 
u G L 1 (R) n L 2 (R), we have 



dfjLu(\) 



Im 



Gg a (x, x', VX)u(x')u(x) 



dx'dxdX 



by Lebesgue's theorem on dominated convergence and polarization yields 

1 



djU«,«(A) 



7T 



Im 



Ge :(7 (x, x' , VA) u(x')v(x) dx'dxdX 



for all »,!)£ L 1 (M) n L 2 (M) since G^ i(T (x, x', VA) is symmetric in x and x' for A > as follows from 
the explicit form (see below). 

2.4. Pointwise decay estimates. As follows from the discussion above, the functional calculus 
for self-adjoint operators yields the representation 
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cos 



TT 



Acos(tA)Im \Ge j( r(x, x', A)l f (x')dx' v{x)dxdX 



for /, v G <S(M) (the Schwartz space), where we have changed variables in the integration with 
respect to A. Our intention is to obtain an expression for [cos (iy'tt^j)/] (x) and thus, we have to 
change the order of integration. However, note carefully that a simple argument based on Fubini's 
theorem does not apply here since the integrals cannot be expected to converge absolutely. In order 
to circumvent this difficulty, first observe that, for any TV G N, we have 



N 



1/N 



Acos(iA)Im \Gg >a (x,x' , X)] f(x')v(x) dx'dxdX < Cn 



which follows immediately from sup x x i eR \G( tCr (x, x', A)| < CV for all A > jj, see Corollary 13.11 
below. Thus, Fubini's theorem yields at least 

o r r r-N 



cos 



— lim 

7T JV->oo 



1/N 



Acos(tA)Im \Gj> >a (x,x' , A)l f(x')dXdx'v(x)dx. 



Next, we distinguish between high and low energies by introducing a smooth cut-off xs satisfying 
Xs(X) = 1 for A G [0, |] and x(A) = for A > 5 where 5 > is sufficiently small. In Sec. [8] below 
we prove the estimate 



(6) 



sup 

x,x'eH 



Acos(tA)Im [G i>a (x,x',X)] (x)- a (x')- a X s(X)dX 



where a £ N and a < 21 + 3. This bound is sufficient to conclude 

"AT 



lim 




lim 

N->oo 



1/N 
N 

1/N 



Acos(tA)Im \ Gi a (x, x', A)] xsWf( x ')dXdx'v(x)dx 



Acos(iA)Im \Gi j0r {x, x', A)l xsWf(x')d\dx'v(x)dx 



by dominated convergence since /,»£ S(R). For the large energy part we show in Sec. [9] that, for 
any a € No, 



(7) 



sup 



lim 

iV-s-oo 



A' 



1/jV 



Acos(iA)Im [G v (x,x',A)] (a?)- a (ar')- a 0(a?')[l - xs(X)}d\dx' 



< 



(ty a / {W{x')\ + W)\)dx' 



which, by dominated convergence, implies 



lim 

N^oo 




lim 



N 

1/N 
N 

1/N 



Acos(tA)Im [Ge :(7 (x, x', A)] [1 — xsW]f(x')dXdx'v(x)dx 



Acos(iA)Im \Gi )(T (x, x', A)l [1 — x&W}f{ x ')dXdx'v{x)dx. 



By adding up the two contributions and using the density of <S(R) in L 2 (R), we arrive at the 
representation 

r-N 

Acos(iA)Im [Gg >a (x, x', A)] d\f(x')dx' 



1/N 



-i 2 

cos(tJ% la )f (x) = lim 

J 7T 7V->oo 

for / G <S(R) and the estimates ©, (J7J) imply the bound 



\w a cos{ty/ne,<T)f\\j 



< 





r 


+ 


/ 










W a 


IACS) 







for 1 < a < 21 + 3 where u> Q (x) := (x) a . An analogous derivation applies to the sine evolution 
and therefore, the proof of our result reduces to oscillatory estimates of the type ([B]) and ©. 

2.5. The main theorem. The main result proved in this work is the following. 

Theorem 2.1. Let (I, a) £ {(0, 0), (0, -3), (1, -3)}, a E N, 1 < a < 2£+3 and set w a {x) := (x)~ a . 
Then the solution operators for the Regge-Wheeler equation satisfy the estimates 



\w a cos(ty/n^)f\\ Lao /m < Ct )Ct {t) a 





f 


+ 


/ 










Wa 






mm) 



and 



sm(tJUiTa) 



< C, a {t) 



-a+1 



for all t > and initial data /, g such that the right-hand sides are finite. 

Remark 2.1. As usual, we prove Theorem 12. II for Schwartz functions f,g G S( 
is then obtained by a standard approximation argument. 
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The general case 



Remark 2.2. For the convenience of the reader we make the behavior of the initial data near 
the horizon more explicit by transforming back to the Schwarzschild r-coordinate. Recall that 
x = r + 2Mlog(2j7 — 1) which implies (x) a ~ |log(2jy — 1)1° as r — > 2M+. Thus, the polynomial 
weights in x translate into logarithmic weights in r. Moreover, we have dx = (1— ^y-)~ 1 dr and this 
shows that the integrability condition near the horizon for initial data g(x) = g{r(x)) transforms 
as 

/o rro ^ r 
\g(x)\(x) a dx ~ / \g(r)\ \log -l)\ a m 
-oo J2M 1 — 

where r = tq corresponds to x = 0. 

We remark that our proof actually applies to more general situations like the analogous problem 
in Hofava-Lifshitz gravity, cf. [33]. The only requirements on the potential are the asymptotics of 
Corollary 12. 1\ the nonexistence of bound states and the nonexistence of a zero energy resonance 
(see Definition 16.11 below) . 



3. Basic properties of the Green's function 

For the convenience of the reader we discuss some well-known properties of the Green's function 
(cf. [15], [56]). 

3.1. The Jost solutions. Recall that the Green's function is constructed with the help of the 
Jost solutions f±(-,z) which are defined by T-in, a f±{-, z) = z 2 f±{-,z) and the asymptotic behavior 
f±(x,z) ~ e ± iZX as x ±00. First we prove that the Jost solutions exist and that they are 
continuous with respect to z in C + \{0} where C + := {z G C : Imz > 0}. 

Lemma 3.1. For every z £ C+\{0} there exist smooth functions f±(-,z) satisfying 

'He,af±(-,z) = z 2 f±(-,z) 

and f±(x,z) ~ e ±lza: for x — > ±00. Furthermore, for every x G R, the functions f±{x,-) and 
/±(av) are continuous in C + \{0}. 

Proof. We only prove the assertion for _/+ since the proof for /_ is completely analogous. The 
variation of constants formula shows that m+(x, z) := e~ lzx f+(x, z), if it exists, satisfies the integral 
equation 

/•oo 

(8) m + (x,z) = l+ K(x,y,z)m + (y,z)dy 

J X 

where K(x, y, z) = (e 2tz< yy- x ) — l) V^ a {y). Conversely, if we can show that Eq. ([8]) has a smooth 
solution, we obtain existence of the Jost solution. However, Eq. ([8]) is a Volterra integral equation 
with a kernel satisfying 

/"OO £j roo 

/ sup \K(x,y,z)\dy < -- , / sup \d x K(x,y, z)\dy < 1 

J a xd{a,y) 1*1 J a x£(a,y) 

for all z G C+\{0} and any fixed a G M (see Corollary 12. ip and thus, Lemma IB. II implies the 
existence of a unique solution m+(-,z) satisfying ||m + (-, 2^) (a,oo) ^ e c ^ z K Furthermore, for 
fixed z G C + \{0}, we have 

/•oo 

/ sup \d%K(x,y,z)\dy < C k 

J a x£(a,y) 
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for all k G No and thus, Lemma lB.21 shows that m + (-,z) is smooth. For the continuity of m + (x, 
fix x G (a, oo), z G C+\{0} and note that 



z + h) - m+(x, z) = g h (x, z) + / K (x, y, z + 7i)[m + (y, 2 + h) - m+(y, z)]dy 

J X 

where 

/•oo 

g h (x,z):= / [K(x,y,z + h) - K(x,y,z)]m+(y,z)dy. 



Now observe that ||<?7j(-,2)||.z/>°(a,oo) -> as /i -> since ||m+(-, 2)||z,oa( 0i00 ) < e c l\ z \ and hence, 
Lemma IB. II implies 

\m + (x,z + /i) - m+(a?,z)| < ||^(-, z)\\ L ^( ai0o) e C/ ^ -> for /i -> 

which shows continuity of m+(ic, ■) in C+\{0} as claimed. For the continuity of m',(x, •) simply 
observe that 

/•oo 

m' + (x,z) = d x K(x,y,z)m + (y,z)dy 

J X 

and the right-hand side of this equation is obviously continuous in z. □ 

3.2. The Wronskian W(f-(-, y/z), /+(•, \/z)). Having established existence of the Jost solutions 
we can now construct the Green's function and the standard procedure yields 

y^Z+Qc, v^eQE - x') + /-(a;, y / i)B(x / - x) 

Gl ltr {x, X , y/Z = — — 

for Imz > where G denotes the Heaviside function. Clearly, W(f-(-, y/z), /+(•, y/z)) / if 
Imz > since otherwise f-(-,y/z) would be an eigenfunction with eigenvalue z contradicting the 
self-adjointness of 7^ CT . However, it is not a priori clear whether the limit Gi a {x, x', y/z) for 
Im.y/z — > 0+ exists. The following observation shows that problems can only occur at z = 0. 

Lemma 3.2. Let A > 0. Then the limit 



W(\/\):= lim T^(/_(-,v / AT^),/ + (-,v / A+^)) 

£-3>0+ 



exists and is nonzero. 



Proof. For brevity we write W(/_, f+){z) instead of W (/_(-, y/z), /+(•, y/z)) and likewise for other 
Wronskians. By Lemma 13. II we know that W (/_,/+) is continuous in C+\{0} and hence, W(y/~X) 
exists for any A > 0. Observe that f'±{x,y/X) ~ iiy/Xe^ 1 ^ 30 for x — > ±oo which follows imme- 
diately from the integral representation in the proof of Lemma 13.11 Thus, W(/+,/+)(A) = 2iy/~X 
which shows that f + (-,y/X) and f+(-,y/~X) are linearly independent for A > 0. Hence, there exist 
A(X) and B(X) such that v 7 ^) = A(A)/+(a?, \/A) + B(X)f + (x,y/X). We conclude 

(9) 2tVX = W(f.,JZ)(X) = W(Af + + BU,AU + Bf + )(X) = -2iVX\A(X)\ 2 + 2sVA|B(A)| 2 
which implies |-B(A)| 2 > 1. However, we have 

(10) W(^X) = W(f-,f+)(X) = W(Af + + BU, f + )(X) = 2iVXB(X) 

and thus, |W(\/A)| > 2y/X which finishes the proof. □ 
Corollary 3.1. The limit 

Gt a (x,x', vA) = lim Ge a (x,x',\/X + ie) 

exists and satisfies 

sup \Ge )0 (x,x , ,VX)\ < C 
li 



for all A > Ao where Xq > is arbitrary. 

Proof. Prom Lemma [3. 2 1 and the asymptotic behavior of the Jost solutions we immediately conclude 



sup 

x'<0,x>0 



W(f-(;V\)J+(;VX)) 



< 1 



for all A > Ao > 0. For the remaining cases use reflection and transmission coefficients ^4(A), B(X) 
(see the proof of Lemma [372]) to express f± in terms of f T and f^. The asymptotic behavior of 
A(X) and B{\) for A — > oo is given by Eqs. (110p and ([9]) and the claim follows. □ 



4. Perturb ative solutions for |xA| small 
In this section we obtain approximations to solutions oirLg^f = A 2 / for |xA| small. The solutions 



0. 



are constructed by perturbation in A around A = 0. We closely follow 

4.1. Zero energy solutions. We first consider zero energy solutions, i.e., solutions of Ti£ }CT f 
By setting v(r) := f(x(r)), the eigenvalue problem rlg^f = A 2 / is equivalent to 

-l 



(11) 



1 



2M\ 
r / 



2M 



-v + 



£(£ + 1) 



2Ma\ 

r 3 J 



A 1 



2M 
r 



and it turns out that for A = this equation can be solved by special functions which will be useful 
later on. However, the following result describes a fundamental system for Hi^f = without 
making use of explicit solutions. 

Lemma 4.1. There exist smooth functions Uj satisfying %e,^u,j = for j = 0, 1 with the bounds 
u {x) = (2£ + l)- 1 x e+1 (l + 0(x- 1+e )) andui(x) = x- e (l + 0(x- 1+e )) for all x > 1 where e G (0,1) 
is arbitrary and the O-terms are of symbol type. The Wronskian is W(uq,ui) = —1. 

Proof. Suppose for the moment that the solution u\ exists and define the function a by u\(x) = 
x (1 + a(x)). Then Hi^ui = is equivalent to 



(12) 



a "f x ) _ — a '(x) 
x 



Vr 



t(l + 1) 



(l + a(x)). 



Viewing this equation as an inhomogeneous equation for a', applying the variation of constants 
formula and integrating by parts, we obtain the integral equation 



(13) 



a(x) 



1 



2^ + 1 



(y 



21+1 



x 2e+1 )y~ 2i 



ViAv) 



£(£ + 1) 



y 



[l + a(y))dy. 



Therefore, if we can show that Eq. (|13h has a smooth solution, we obtain existence of u\. However, 
Eq. (|13p is a Volterra integral equation of the form 

/■oo 

K(x,y)dy+ / K(x,y)a(y)dy 



a(x) = 

with a kernel K satisfying \d^d l y K(x, y)\ < Ckjy~ 2+£ ~ k ~ l for 1 < x < y and an arbitrary e G (0, 1) 
(cf. Corollary 12 . X [) . Therefore, Lemma IB.2I implies the existence of a unique smooth solution 
a G L°°(l,oo) and Eq. (|13p shows that in fact |a(x)| < x~ 1+£ for x > 1. Furthermore, the first 
derivative of a is given by 



a'(x) 



d x K(x,y)(l + a(y))dy 



and this implies |a'(x)| < x 2+£ for x > 1. The estimates for the higher derivatives follow from 
Eq. (I12p . the Leibniz rule and a simple induction. 
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For the second solution uq we use the Wronskian condition 

l = W( Ul ,u )= (^Y«i 

which yields 

u Q {x) = ux(x) / u{ 2 {y)dy 



j-i 



where x\ > is chosen so large that u\(x) > for all x > x\. Note that u^ 2 (x) = x 2e (l + b(x)) for 
a suitable b satisfying |&( fc )(x)| < Ck{x)~ 1+£ ~ k (apply the Leibniz rule and Lemma I A. 1|) . Inserting 
the asymptotic expansion for u\ yields uq{x) = (2£ + + 0(x~ 1+£ )) where the O-term 

behaves like a symbol. □ 

We construct another pair vo, v\ of zero energy solutions with specific asymptotic behavior as 
x — > — oo. This is considerably easier than the above construction for the solutions Uj due to the 
exponential decay of the Regge- Wheeler potential V^ CT (x) as x — > — oo. 

Lemma 4.2. There exist smooth functions Vj for j = 0,1 satisfying Tin^Vj = and v$(x) = 
x{l + 0(x^ 1 )) as well as v\(x) = 1 + 0(x~ 1 ) for all x < — 1 where the O-terms behave like symbols 
under differentiation. 

Proof. For x < — 1 consider the Volterra equations 



and 

v\{x) = l-l (y - x)Vp a {y)vx(y)dy 



which have smooth solutions according to Lemma lB.21 since the potential Vi tCT (y) decays exponen- 
tially as y — > — oo (see Corollary 12. ip . Obviously, we have the asymptotic behavior vo(x) ~ x, 
vi(x) ~ 1 as x —7- — oo, Vj behave like symbols under differentiation and satisfy T-Le )CT v j = 0, as a 
straightforward calculation shows. □ 

Corollary 4.1. The solutions uj for j = 0,1 can be uniquely extended to all of M and we have 
Uj(x) = 0(x) as x —> —oo where the O-term is of symbol type. 

Proof. Since the potential Vi tCT is smooth on R, the solutions Uj, originally defined on [l,oo) only, 
can be uniquely extended to all of M. by solving appropriate initial value problems. Since the 
solution pair vq, v\ forms a fundamental system for the equation Hi^f = 0, Uj can be written as a 
linear combination of vq, v\ on (— oo, —1] and everything follows from Lemma 14.21 □ 

4.2. Construction of the perturbative solutions. Next, by perturbing in A around A = 0, we 
obtain useful approximations to solutions of Hi tCT f = A 2 / for \x\\ small. 

Lemma 4.3. There exist constants xo,Xo,5 > and smooth functions Uj(-,X) satisfying 

He >a uj(-,\) = X 2 Uj(-, A) 

for j = 0,1 and W(uo(-, X),u%(-, A)) = — 1 such that Uj(x, A) = Uj(x)(l + aj(x, A)) where 

\a,(x,X)\<x 2 X 2 

for all X € (0, Ao) and x E [xc^A^ 1 ] provided that £ > 1. In the case I = we have the weaker 
bounds 

|ao(x,A)| < x 2 A 2 and \a\(x,X)\ < xX 

in the above ranges of x and X. 
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Proof. Let xo > be so large that uq(x) > for all x > xo- A straightforward calculation shows 
that, if the function h solves the integral equation 



(14) 



h(x,X) = 1 - A 2 



uo{y)ui(y) - ul{y) Ul ' 



u (x) 



h(y,\)dy, 



then uq(x,\) := uo(x)h(x, A) satisfies He^uo^, A) = A 2 uo(-,A). Eq. (fT4"|) is a Volterra integral 
equation 



/i(x,A) = 1+ / K{x,y,X)h(y,X)dy 

where the kernel is of the form 

K(x, y, A) = A 2 \y(l + 0{y- 1+£ )) - x^ y 2e+ ' ! (1 + 0(*- 1+£ ) + 0(y- 1+£ )) 



for x,y > xq and the O-terms are of symbol type (see Lemma l4.ip . This shows \d^d l y K(x, y, A)| < 
Ck iy 1 ~ k ~ l \ 2 for xq < y < x and hence, 



/ sup |^(a!,y,A)|dy<l 

Jxq {x:xo<y<x} 



for A € (0, Ao) and all k G No where Ao > is arbitrary. Thus, Lemma lB.21 shows that Eq. (|14p has 
a unique smooth solution h(-,X) satisfying \\h(-, X)\\Loor XOt x-i-\ < C for all A G (0, Ao). Therefore, 
Eq. flUD implies \h(x, A) - 1| < x 2 A 2 for A G (0, A ) and x < x < A -1 . 

Now choose <5 > and Ao > so small that \h(x, A) — 1| < ^ for all A G (0, Ao) and x G [xo, <5A -1 ]. 
We use the Wronskian condition 



-l = W(uo(;X), Ul (;X)) 

to construct the second solution u±(x, A), i.e., 



A) V 
u (-, A) 



u 2 (,A) 



«i(x, A) = uq(x, A) 



5A" 



u 2 {y,X)dy 



which implies 



(15) 



ai(x, A) 



uq(x) 
ui(x) 



(l + a (x,A)) / u '(y)(l + a (y,X))dy -1 



where ao(y, A) := (1 + a^{y, A)) 2 — 1 inherits the bound of ao- Now, by inserting the asymptotics 
of uq from Lemma 14.11 we obtain 



for x — > oo and this shows 



u o( x ) / u (y)dy ~ x~ 



ni(x) = u Q (x) / n (y)c?y 

•J 
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since u± is uniquely determined by the asymptotic behavior u\(x) ~ x e for x — > oo (cf. Lemma 
I4.ip . Using this, Eq. (fT5|) and the asymptotics of ao, cto> we obtain 



ai(x,X) 



+ 



up(x) 
ui(x) 

u (x) 



u 2 (y)dy + a (x,X) 



sx- 



u 2 (y) d y 



1 + a {x, A)) 



s\- 



ui{x) 
= 0(x 2e+1 
= 0(x 2 X 2 ) + 0((xX) 2e+1 ) 
which implies the claim. 



u 2 (y)a (y,\)dy 



Q{X 2l+l ) + 0{x 2 X 2 ){0(X 2l+1 ) + 0{x~ 2l ~ 1 )) + 0(x- 2£ - 1 )0(x 2 X 2 ) 



□ 



4.3. Estimates on the derivatives. Next, we study derivatives of the above constructed solu- 
tions. 

Proposition 4.1. The functions aj for j = 0,1 from Lemma \4.3\ are of symbol type, i.e., 

\d^ aj (x,X)\<C k , m x 2 - k X 2 - m 

for k,m G No and X G (0, Ao), x G [xo,<5A -1 ] provided that £>1 where xq,Xq,5 > are constants. 
In the case £ = we have the weaker bounds 

\8£8?ao(x,\)\ < C Km x 2 ~ k X 2 ~ m and \d k x d^ ai {x,X)\ < C^x 1 ^ X 1 '™ 

in the above ranges of x and X. 

Proof. We use the notations from the proof of Lemma 14.31 and proceed by induction. We have to 
consider mixed derivatives and therefore, we order the set No x No according to 

(0,0), (1,0), (0,1), (2,0), (1,1), (0,2),... 

which defines a bijection n : N x N -)■ N by n(0,0) = 0,n(l,0) = l,n(0, 1) = 2, etc. Now 
fix (k,m) G N x N and assume that |c^,<9^(/i(x, A) - 1)| < Ci t jX 2 ~ l X 2 ~i holds for all with 
n(l,j) <n(k,m) and A G (0, Ao), x G [xo,A -1 ]. We have to show that this implies \d k d^ 1 h(x, A)| < 
Ck',m' x2 ~ k X 2 ~ m , where n(k',m') = n{k,m) + 1. There are two possibilities: Either (k',m') = 
(m + 1,0) (if k = 0) or (k',m') = (k — l,m + 1). In the former case we have with Ki{x, A) := 
d l x K(x,y,X)\ y=x = Oix^X 2 ) (cf. Lemma El, 

d™ +1 h(x, X) = jr d™- 1 [mix, X)h(x, A)] + f d™ +l K(x, y, X)h(y, X)dy. 



1=0 



Now observe that by assumption \d x n 1 [ki(x, X)h(x, A)] | < x 2 ( m+1 )A 2 for A G (0, Ao) and x G 
[xq, A -1 ] and hence, |9™ +1 /i(x, A)| < x 2 ~( m+1 )A 2 . In the latter case we have, provided k > 2, 



d™ +1 d h x - l h(x,X)\ 



^2c^ +1 dt 2 - l Mx,X)h(x,X)}+ I* d™ +1 k fc - 1 ^(x,y,A)%,A) 

J_n Jxr, 



1=0 Jx ° 
< x 2-(fc-l) A 2-(m+l) 
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dy 



by assumption. If k = 1 we have 

min{m+l,2} 



d™ +1 h{x,X) = £ ) diK(x,y,\)d™ +1 ->h(y,\)dy 

j = l \ J / Jxq 

+ [* K(x,y,X)d™ +1 h(y,X)dy 

JXQ 

and thus, by assumption, the derivative d™ +1 h(x, A) satisfies a Volterra equation of the form 
d? +1 h(x,\) = 0(x 2 X 2 ~^) + [ X K(x,y,\)d™ +1 h(y,\)dy 

and the basic estimate from Lemma IBTTI yields \d™ +1 h(x, A)| < x 2 A 2 ~^ m+1 - ) for all A G (0,Ao) and 
x G [x , A -1 ]. 

For the second solution we use the representation 

Uq ( x ) ( f°° rSX- 1 \ 

(16) ai(x,A) = — -- - / UQ 2 (y)dy + a (x,\) u 2 {y)dy\ 



ui(x) 



■uo(x) /' 5A_ 
H -—-(l + a (x,X)) UQ 2 (y)d (y,X)dy 



ui(x) 

from the proof of Lemma 14.31 where, as before, do(x, A) := (1 + ao(x, A))~ 2 — 1. Lemma lA.ll and 
the Leibniz rule show that do inherits the bounds of o,q, i.e., \dxd™ao(x, A) | < Cfc )m x 2_ ' c A 2_m for 
all fc,m E No and A G (0, Ao), 2; G [xo, A -1 ]. Thus, all functions on the right-hand side of Eq. (fTB]) 
behave like symbols under differentiation with respect to x and A. Therefore, as in the proof of 
Lemma l4.3( we have 

ai {x, A) = 0(x 2 A 2 ) + 0{{x\f t+l ) 
where the O-terms are of symbol type which finishes the proof. □ 

4.4. Refined bounds for A— derivatives. As a next step we prove a refinement of the estimates 
for the solution no which shows that we can effectively trade A -1 for x in the bounds for the 
A-derivatives of ao- 

Lemma 4.4. The function ao(-, A), defined by uo(x, A) = uo(x)(l + ao(x, A)), satisfies the estimates 

|^ m a (x,A)| < C m x 2m and |d 2m+1 a (x, A)| < C m x 2m+2 A 

for all A G [0, Ao], x G [xo, A -1 ] and m G No where Ao > is a sufficiently small constant. 

Proof. We use the notations from the proof of Proposition 14.11 and proceed by induction. The case 
m = has already been proved in Proposition 14.11 Now fix m G N and assume that \d 2 x j (h(x,X) - 
1)1 < CjX 2 i and \d 2j+1 (h(x, X)-l)\ < CjX^ +2 XioT j < mm the above range of A and x. According 
to the proof of Proposition 14.11 we have 

d 2m h(x,X) = ( 2 ™ ) f d 2 K(x,y,X)d 2 x m - 2 h(y,X)dy + 2m f d x K(x,y, X)dl' n - 1 h(y, X)dy 
\ / Jxo J XO 

+ [ X K(x,y,X)d 2 x m h(y,X)dy 
Jx 

= f 0(y)0(y 2m - 2 )dy + f ' 0(yX)0(y 2m X)dy + f K(x,y, X)d 2m h(y, X)dy 

J Xo J XO J Xo 

= 0(x 2m )+ [ X K(x,y,X)d 2m h(y,X)dy 

J XQ 
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by assumption and the estimate in Lemma [B.ll implies \d 2m h(x, X)\ < x 2m for A G [0, Ao] and 
x E [xq, A -1 ]. For the odd derivatives we proceed analogously and obtain 



d 2 x m+1 h(x,X) = ( 2m 2 +1 ) jyiK(x,y,X)dl m - 1 h(y,X)dy 



+ (2m + l) d x K(x,y,X)d 2m h(y,X)dy + K(x,y, X)d 2 x m+1 h(y, X)dy 

0(y)0(y 2m X)dy + f 0{yX)0{y 2m )dy + [* K(x,y,X)d 2 x m+1 h(y,X)dy 

x Jx Jx 



0(x 2m+2 X) + [ X K(x,y,X)d 2 x m+1 h(y,X)dy 

JXQ 



by assumption and again, Lemma IB . 1 1 yields the claim. □ 

Note that, by construction, we have uo(xq,X) = uo(xo) and u' (xq,X) = u' (xq) (see the proof 
of Lemma |4.3|) which shows in particular that uo(xq, A) and u' (xo,X) are smooth functions of A. 
Next, we prove similar bounds for the function a\ but unfortunately, the situation here is a bit 
more complicated. 



Lemma 4.5. Let ui(x,X) = ui(x){1 + a\(x, A)) and xq,Xq,5 > be as in Lemma \4-3[ Then, for 
all X G (0, Ao) and all x G [xo,<5A _1 ], we have the estimates 

\d 2 x m ai (x,X)\ < C m x 2m and |d 2m+1 aiO, A)| < C m x 2m+2 X 

provided that m < £ — 1. Furthermore, for higher derivatives we have the bounds 

\d 2 x e+m ai (x,X)\<C m x 2e X- m 

for m G No in the above ranges of x and X. 

Proof. The function a\(x, A) is given by 

ai(x, A) = ^44(1 + a (x, A)) u Q 2 (y)(l + 5 (y, X))dy - 1 



h(x, A) := r «o 2 (y)(! + 5o(v» " L 



with ao(x, A) = (1 + ao(x, A))~ 2 — 1, see Eq. (fT5]> . Thus, in view of Lemma l4~4l it suffices to prove 
the claimed bounds for 

m(x) 

For k G N we have 

l 

d\h(x, A) = (^(^^-^'(l + ao^A))!^^ 2 ) 

i=o 

+ ^ T u 2 (y)d k x (l + d (y,X))dy 



ui(x) 

and, by using the symbol behavior of ao, we infer 

E ^ («o "(tfA- 1 )^" 1 -^! + ao(y, A))^-.^ 2 ) = 0(* 2 ^)0(A 2 ^- fc ). 

If A; = 2m we have 0(x 2m )0(A 2m_A: ) = 0(x 2m ) provided that m < I. If /c = 2m + 1 we obtain 
0(x 2m )0(A 2m - fc ) = 0(x 2e+1 - 2e+2m+1 X) = 0{x 2m+2 X) 
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provided that m < £ — 1. In order to estimate the integral term note that ao inherits the bounds 
of ao from Lemma 14.41 as a consequence of the Leibniz rule and Lemma IA.1I For even derivatives 
we therefore have 

uo(x) ^ ^ 



ui(x) 



we use d 2 x m+1 d (y,X) = 0(y 2m+2 X) and infer 



u \y)dl m (l + ao(y, X))dy = 0(x 2e+1 ) / 0{y~ n ~ 2+2m )dy 

J X 

= 0(x 2£+1 )(0(\ 2i+1 - 2m ) + 0(x- 2e - 1+2m )) 
and, provided that m < £, we obtain 0(x 2e+1 )0(X 2iJrl ~ 2m ) = 0(x 2m ) as before. For odd derivatives 

7 A 
Uq{x) 



ui(x) 



u» 2 (y)d 2 x m+1 (l + a (y,\))dy = 0(x 2 ^) / 0(y- 2 ^ 2m+2 X)dy 

J X 

= 0{x 2t+1 ){0{X 2t - 2m ) + 0(x- 2e - 1+2m+2 X)) 



and 0(x 2£+1 )0(X 2e 2m ) = 0{x 2m+2 X) provided that m< £—1. The claim for the higher derivatives 
follows directly from the symbol behavior of the above O-terms. □ 

At this point it is convenient to introduce a new notation. 

Definition 4.1. For N G No we write fix) = 02n(X) if, for a constant a > 0, 

(1) / : (0,o) -> R is smooth, 

(2) |/ (fc) (^)l < Cfc for k<2N and all x G (0, a), 

(3) |/ (2JV+fe) (z)| < C k x~ k for all fc G N and x G (0,o), 

(4) lim a; _ > o+ / (2 * -1) (a0 = for all 1 < A; < N. 
Similarly, we write f(x) = O^N+lix) if, for a constant a > 0, 

(1) / : (0,a) -»-R is smooth, 

(2) \f {k \x)\ < Cfc for k < 2N + 1 and all x G (0,a), 

(3) |/( 2Ar+1+fe )(x)| < C k x- k for all A; G N and x G (0,a), 

(4) lim x ^ + / (2fc) (>) = for all < k < N. 

Note carefully that the crucial difference between O2N and 02jV+i is in condition (4). We also use 
the symbols O2N+1 and O2N to denote generic real-valued functions with the respective properties. 



Corollary 4.2. Let u\(x, A) and xq, Xq > be as in Lemma \4^[ Then we have ui(xq, A) = O2HI) 
and u'i[xq, X) = 02^(1) for all X G (0, Ao). 

Proof. The first assertion ui(xq,X) = 02^(1) follows immediately from Lemma 14.51 For the second 
one note that Uo(xq) / and by construction we have —1 = Uo(xo)u' 1 (xo, A) — u' (xq)u(xo, A) which 
implies 

u^xq, A) = -y — r = O2H 1 ) 

uo {xq) 

by the first part. □ 

4.5. Extension of v,j(x, A) to negative values of x. Finally, we extend the solutions Uj(x,X) to 
negative values of x and prove appropriate estimates. 

Lemma 4.6. The functions Uj(x,X), j = 0,1, from Lemma \4-3\ can be smoothly extended to 
x G [— A , xq] for X G (0, Ao) where Ao > is a constant. Furthermore, in the above ranges of x 
and X, the function uq satisfies the bounds 



\d^u (x,X)\<C m (x 
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m+l 



m G No, whereas for u± we have the estimates 

|^«i(x,A)| <C m (x) m+1 

if m < 2£ and 

\dl i+m Ul (x,\)\<C m (x) 2l+1 \- m 

for m G No- 

Proof. The coefficients of the equation 7i a ,if = A 2 / are smooth on R and thus, any solution of 
T~L a ,tf = A 2 / can be smoothly extended to all of R by solving an appropriate initial value problem. 
Applying the variation of constants formula and noting that uq(xo, A) = uo(xq), u' (xq, A) = u' (xo) 
shows that the solution uo(-, A) satisfies the integral equation 

f-XQ 

uq(x, A) = u (x) + A 2 / [v (x)vi(y) -v (y)vi(x)]u (y,\)dy, 

J X 

see also Corollary I4.lt where Vq, x>\ are the smooth extensions to (—00, xq] of the functions con- 
structed in Lemma [4.2i This is a Volterra equation with a kernel K(x,y,X) := \ 2 (vo(x)vi(y) — 
vo(y)vi(x)). According to Lemma l4~2l we have the bound \K(x,y, A)| < A 2 ((x) + (y)) which implies 

/XQ 
sup \K(x,y,X)\dy < 1 
-A- 1 a:6(-A- l ^so) 

and hence, Lemma [B.ll and Corollary 14.11 show that |uo(x,A)| < (x) for all A G (0, Aq) and x G 
[— A _1 ,xo] where Ao > is the constant from Lemma 14.31 We proceed by induction. Fix m G No 
and assume that we have \d l x uo(x, A)| < Ci{x) l+1 for alH < m in the above range of A and x. This 
implies 

min{m+l,2} , , 

d™ +1 u (x, A) = ^ \ l )l y ' A ) 5 r +1_ ^o(y, A)dy 

+ / K(x,y,\)d™ +1 u (y,\)dy 

J X 

min{m+l,2} 



2 Q(A 2 - ; (x) 2 ) 0((x) 1+m + 1 -') + J K(x,y,\)dT +1 u (y,\)dy 



1 1 o{(xY) 

and the estimate from Lemma IB . 1 1 yields \d™ +1 uo(x, A)| < (x) m+2 as claimed. Note carefully that 
the index I in the last sum is at most equal to 2 and therefore, we only estimate nonnegative powers 
of A here. 

For the second assertion we proceed similarly and note that a straightforward calculation as well 
as the variation of constants formula show that u\(-, A) satisfies the Volterra equation 

ui(x ,X)v' l (x )-u , l (xo,X)vi(x ) ui(x , X)vq(xo) - u'iixo, X)vo(x ) 

Ul{X,X) = — r V {X) — r Vl(x) 

W(vo,vi) W(v , Vl ) 
+ / K(x,y,X)ui(y,X)dy. 



According to Corollary 14.21 we have 

ui(x,X) = O 2 £(l)v (x) +0 2 e{l)vi(x) + / K(x,y,X)m(y,X)dy 

J X 

and, since Vj(x) = 0((x)), we obtain |tti(», A)| < (x) for all x G [— A -1 ,^] an d A G (0, Ao) by 
Lemma iB.li Again, we proceed by induction and first we consider the case m < 21. If £ = there 
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is nothing left to prove, so we restrict ourselves to t > 1. Assuming that \d l yU\{x, A)| < (x) l+1 holds 
for all Z < m and a fixed m < 21 — 1 , we infer 

min{m+l,2} 

a™ +1 Ul (x,A) = 0((x))+ J] { I )J d[K(x,y,X)d^ +1 - l u 1 (y,X) dy 



i=i 

■I'D 



0((a:»0((tf>"») 



+ / K(x,y,\)d™ +1 Ul (y,\)dy 



= 0«*> m+2 ) + f ° K(x, y, \)d™ +1 ui(y, X)dy 

J X 

and Lemma lETTl yields \d™ +1 ui{x } A)| < (x) m+2 for all x G [A _1 ,a;o] and A G (0, A ). The claim for 
the higher derivatives follows by a similar induction. □ 

5. Perturbative solutions for |xA| large 

The solutions u,(-,A) obtained by perturbing in energy cannot directly be matched with the 
Jost solution / + (-,A) since the approximations in Lemma 14,31 are valid for \x\\ small whereas the 
behavior of f+(x, A) is known only for fixed A and x — > oo. Thus, we construct another set of 
solutions to 7ii. a f = A 2 / by perturbing the potential. 

5.1. Construction of the perturbative solutions. To do so, we first rescale the equation 
Htcrf = ^ 2 f by introducing a new independent variable z := Xx. Setting f(z) := f(X~ 1 z), 
the equation Hi^f = A 2 / is equivalent to 

(it) /" + (i - e -^±^j / = x-^a^z)! 

where Ui a {x) := Ve a (x) — ^+D. Now observe that the equation 

r + (i-^)/=. 

has the fundamental system {y/zJ£+u2(z), ^/zY^^z)} where Jg+i/2 and are the Bessel 

functions pQ. Thus, for a small right-hand side of Eq. (|17p we expect to obtain solutions of Eq. (|17p 
by perturbing the Bessel functions. According to Corollary 12.11 the right-hand side satisfies the 
estimate X~ 2 Uz^{X~ l z) < A _2 (A~ 1 2:) _3+£ = A~ 2 x _3+£ for x — > oo. Thus, our approximation is 
expected to be good if A~ 2 x -3+e is small. Smallness can be achieved by fixing A > and letting 
x — > oo which is required for the matching with the Jost solution /+. However, we can also enforce 
smallness by setting x = A~ 1+e (for a small e > 0) and letting A — > 0. For A > sufficiently small, 
we have |xA| < 5 and the matching with the solutions Uj(-, A) can be done as well. As a result, we 
obtain a good approximation to the Jost solution f+(x, A) at a finite x for A — > 0. 

Lemma 5.1. There exists a smooth solution 4>e(-,X) of Eq. fli7| ) such that 

Mz, A) = (3 e ^H+ +1/2 (z)(l + b t (z, A)) 

where H^ +1 / 2 := -h+i/2 + ^+1/2 ^ s the Hankel function and := zy^e^"2. For all X G (0, 1), the 
function be satisfies the bounds 

MzA^^cx 1 -^ 2 ^ 

for all z G [A £ , 1] and 

Mz,x)\ kcz-^x 1 - 6 

for all z G [1, oo) where e G ^0, *s arbitrary. 
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Proof. Set ipe(z) := Pe^/zH^-^i^z) and observe that |^(#)| > for all z > 0. Furthermore, we 

have W(ipe,lpe) = —2i which follows by noting that ifi(z) ~ e lz and ip' e (z) ~ ie' LZ for z — > oo. A 
straightforward calculation shows that, if h satisfies the integral equation 



(18) 



h(z,X) = 1 



1 

2i 



<pt(y)<pt(v) - fe(y) 



\- 2 U^(X- 1 y)h(y,X)dy, 



then 4>e(z, A) := ipe(z)h(z, A) is a solution to Eq. (fTTj) . Eq. ([18]) is of the form 

poo 

h(z,X) = 1+ / K(z,y,X)h{y,X)dy. 



Recall the asymptotic behavior ipe(z) = ci@£Z £+1 (l + O(z)) + ic2$iz~ IL {\ + 0{z)) for z — > where 
ci,C2 are nonzero real constants and the O-terms are smooth (cf. pQ). Furthermore, we have 
(pi(z) ~ e tz as z — > oo and \X~ 2 U(,^(X~ x z)\ < A 1_e 2: _3+e for all X,z > with, say, A _1 z > 1 (see 
Corollary 12. ip . This shows that, for k G No, 

f — x (fc) 



< a 



for all z > and hence, \d k K(z,y,X)\ < X 1 ~ e y~ 3+e (l + y~ n ) for all A G (0,1) and all y,z with 
X s < z < y. Thus, we have 



/i(A) := 



sup |K(^,i/,A)|dy < A 1 -^ 3 ) 



for all A G (0, 1) and hence, /i := sup Ag ( 0; i) /-*(A) < oo provided that e < ^pg- Applying LemmaEJ 
we conclude that Eq. (fT8l) has a unique smooth solution h(-,X) satisfying \\h(-, A)||//x>r\e )0O ) < e M 
for all A G (0, 1). Thus, Eq. CEHD implies 



|fc(*,A)-l|< / \K(z,y,X)\dy< / A 1-e y -3+e d?/ < A 1-e 2 _2+£ 
for all A G (0, 1) and all z G [1, oo). Similarly, for z G [A e , 1], we have 



\h{z,X)-l\ 



X l-e y -3+s( 1 + y -M )dy 



\K(z,y,X)\dy< 
< X 1 - 6 (V 2+e + z- 2 - 2e+: ) < a j -- i2 '-' !i 

as claimed. 

5.2. Estimates for the derivatives. 



-2l\ 



□ 



Lemma 5.2. For a// A G (0, 1) and m,k G No, f/ie function bi from Lemma I5.il satisfies the 
estimates 

\d k d?b e (z,X)\ < C Km z- k X l ~ e ^- m 

for all z G [A £ , 1] and 

\d k z d?h{z,X)\ < C Km z- 2+£ - k X l - £ - m 
for all z G [1, oo) where e G ^0, 2^3) *s arbitrary. 
Proof. The function h := 1 + b# satisfies the equation 



roo 

h(z,X) = 1+ / K(z,y,X)h(y,X)dy 

J z 
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where 



K(z,y,X) :=4 



\- 2 U t , a {X~ l y) 



and ife(z) := f5iyfzH^ +1 j 2 (z) : see the proof of Lemma I5TT1 Recall the asymptotic behavior <pi(z) = 

e lz {\ + Oc(z -1 )) for z — > oo where the Op-term is of symbol type. On the other hand, we 
have (pi(z) = ci(3£z e+1 (l + 0{z)) + ic2@£Z~ l {l + O(z)) for z — > where c\,C2 are nonzero real 
constants and the O-terms are smooth (see PQ). Thus, by the Leibniz rule and Lemma fA.lt we have 
K(z,y,X) = (1 - e 2i ^- z ^)O c {y- ?J+£ X l - £ ) for all 1 < z < y and K(z,y,X) = Oc(y~ 3 ~ 2e+e A 1 " 6 ) = 
Oc(y~ 1 A 1 ~ £ ^ +3 )) for A e < z < y < 1 where all O-terms are of symbol type. Let n : No x No — > No 
denote the bijection from the proof of Proposition 14.11 As before, we proceed by induction. Fix 
(k,m) G N x N and suppose we have \d l z d{(h(z , X) - 1)| < Ci tj z^ 1 X 1 ^ £< ^ 2e+3 ^ j for all with 
< n(k,m) and A G (0,1), z£ [A £ ,l]- We have to show that this implies |c^?' 9™' (h(z, A) — 1)| < 
Oj./ !m /2 _fe 'A 1_e ^ +3 ^ _m ' where n(k',m') = n(k,m) + 1. If k = we have (k' ,m') = (m + 1,0) and, 



with «j (z, A) := (T Z K (z, y, A) | 



y=zi 
m 



ar +1 ^(z, A) = - J] S^-' [«,(*, X)h(z, A)] + / d? +1 K(z, y, X)h(y, X)dy, 
i=o ^ z 

see LemmaEl By assumption we have \d™~ 1 [ki(z, X)h(z, A)] | < z -(m+i) for A G (0,1), 
z G [A £ ,l] and this implies \&? +1 h(z, A)| < z -(m+l) A l- £ (2^+3)_ If > 1 we have (k',m') = (k - 
1, m + 1) and 

(19) aj^flf-^foA) =-X)ST +1 ^" 2 " , [«i(^A)/i(^A)]+ / [dt'Kiz^AMy^)] dy. 

i=o Jz 

If k > 2, Eq. (HID and the assumption shows that \d%+ 1 6%- 1 h(z, A)| < «-(*-!) A 1-6 ^ 3 )-^ 1 ) for 
A G (0, 1), z G [A e , 1]. If Jfe = 1, Eq. (PSD is of the form 

d™ +1 h(z,X) = o c (A 1 - £ ^ +3 )-( m+1 )) + / K(z,y, X)d™ +1 h(y, X)dy 



by assumption and therefore, Lemma fB . 1 1 yields \d^ +1 h(z, A)| < \ 1 - £ ( 2e + 3 )-( m + 1 ) for A G (0, 1) and 
z G [A £ , 1]. This proves the first estimate. 

For the second estimate we proceed by an analogous induction and write 



o 



K(z,y,X)h(y,X)dy = / K(z,rj + z, X)h(rj + z, A)dr? 

(1 - e 2i ")O c ((»7 + ^)- 3+£ A 1 - £ )^(r ? + z, A)dr? 
for z > 1 where the Oc-term is of symbol type. Thus, the assumption yields 

POD 

d™ +1 h(z,X) = O c (z- 2+£ -( m+1 )A 1 - £ ) + / K(z,y,X)d™ +1 h(y,X)dy 



o 



and Lemma TB.1I implies \d™ +1 h(z, A)| < z 2+e ( m+1 )A 1 6 for A G (0,1) and z > 1. Analogously, 
we obtain 

c^+^/^A) = O c (z- 2+£ -(^ 1 )A 1 - £ -( m+1 )) + / if (z, y, X)d7 b+1 dy~ 1 h(y, X)dy 



and again, Lemma [Bj] yields |Sj _1 ^ +1 /i(z, A)| < z - 2 +e-(fc-i) A l-e-(m+i) for A G ^ 1) and z > 1 
which finishes the proof. □ 
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6. Matching with the Jost solutions 

In this section we match the Jost functions /±(-,A) with the solutions Uj(-,X) which allows us 
to calculate the asymptotic behavior of /±(-,A) and W(f- (•, A), /+(•, A)) in the limit A — > 0+. 

6.1. Matching with /+(•, A). Note that the solution <p£ constructed in Lemma 15.11 has the as- 
ymptotic behavior <j)g(z,X) ~ e lz for z — > oo which shows that f + (x,X) = 4>i(Xx,X). Thus, we 
have found a representation of the Jost solution A) which is valid for all A E (0, 1) and all x 

with A~ 1+£ < x < oo with a sufficiently small e > 0. For given e, 5 > we can always accomplish 
A~ 1+e < (5A _1 for all A E (0, Ao) provided Ao > is chosen small enough. Thus, at x = A _1+e for A 
sufficiently close to 0, we can match the solutions / + (-,A) and Uj(-,X) (see Lemma l4T3|h 



Lemma 6.1. The Wronskians Cj(X) := W(f+(-,X),Uj(-,X)) for j = 0,1 have the asymptotic 
behavior 

c+(A) = iaoPtX-tQ. + 0(A e ) + iO(X £ ^ +2 ^)) 

and 

C +(A) = ai fc\ e+1 (l + 0(A £ ) + iO(A~ 2& )) 

as A — > 0+ for a sufficiently small e > where ay are real nonzero constants and all O-terms are 
of symbol type. 

Proof. We have f+(x, A) = fyyXxH^ l ^ 2 (Xx)(l + bi(Xx, A)) by Lemma loTTl Note that d x be(Xx, A) = 

Xd z bi(Xx, A) as well as d\bi(Xx, A) = xd z bi(Xx, A) + d\bi(Xx, A) and hence, by Lemma [5\2l and the 
chain rule, we infer 

\d*d^b e (Xx,X)\ < C kim x- k \ 1 -< 2e+s >- m 

for all k,m E No- By the same reasoning we obtain the symbol character (with respect to differen- 
tiation in x and A) of the O-terms in 

faVXx~H+ +1/2 (Xx) = aift(Ax) m (l + 0{Xxj) + ia p e {Xx)- e (l + O(Xx)) 

where ao, ot\ are nonzero real constants. This shows that 

f+(x, A) = (a 1 p e (\x)'+ 1 + iaoM^r") (1 + O c (A 1 - £ ^+ 3 )))(l + O(Xx)) 

where the O-terms are of symbol type and the representation is valid for all A E (0, 1) and all 
x E [A +E ,A ]. Differentiating this expression with respect to x and using the symbol character 
of the O-terms we obtain 

f' + (x, A) = ((£ + l)ai^A(Ax)^ - ifeoAA(Ax)-^ 1 ) (1 + O c (A 1 - £ ( 2£+3 )))(l + 0{Xx)). 
Evaluation at x = A~ 1+e yields 

i+(A- 1+£ , A) = (anPtX'V+V + uxoPtX-**} (l + 0(A e ) + iOCA 1 " 6 ^ 3 ))) 

and 

f' + {X' l+£ , X) = ((£ + l)a 1 p e \ 1+Bt - ilaoPeX 1 -^) (l + 0(A £ ) + ^(A 1 ^ 2 ^)) 
for e > sufficiently small. Furthermore, by Lemma 14.31 we have 

uo(x, A) = u (x)(l + 0(x 2 A 2 )) , u' (x, A) = u' (x)(l + 0(x 2 A 2 )) 

and 

ut(x,X) =u 1 (x)(l + 0(x\)) , u[(x,X) = ui(x)(l + 0(xA)) 
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for all A G (0, Aq) and all x G [xo,<JA _1 ] where \q,5 > are sufficiently small and xq > is 
sufficiently large. Thus, by choosing A sufficiently close to 0, we obtain A _1+£ G [xo,(5A _1 ] an d we 
can evaluate the above expressions at x = A~ 1+e which yields 

u (X~ 1+£ ,X) = (21 + i)- 1 A- £ - 1+e ( £+1 )(l + 0(A 2e )), 

u' (\- 1+e , X) = (£ + 1)(2£ + lr^+^l + 0(A 2£ )), 

Ul (A- 1+e ,A) = A £ - e€ (l + 0(A e )), 

ui(A- 1+e ,A) = -£\ m " e(m) (l + 0(A e )) 

by Lemma 14.11 and all O-terms are of symbol type. The claim now follows from a straightforward 
computation. □ 

6.2. The Jost solution /_(-,A) in the limit A — > 0+. The Jost solution /_(-,A) satisfies the 
Volterra integral equation 

/-(x, A) = e-^' - sin{X{y x - x)) v e Ay)f-(y, m 

as can be seen from the definition and the variation of constants formula. The decay properties of 
the potential Vg. CT are crucial for the behavior of f±. Since V£. CT decays exponentially as x — > — oo 
(Corollary 12. ip . the situation for /_ is much simpler. In fact, /_ behaves essentially as in the free 
case Vt )tr = 0. 

Lemma 6.2. Let a G R and Ao > 0. Then the Jost solution f-(x, A) = e~ tXx m-(x, A) exists for all 
A € [— Ao, Ao] and, for m G No, we have the bounds \d™m-(x, A)| < C m as well as \d™m!_(x, A)| < 
Cm for all x G (— oo,a] and a// A G [— Ao, Ao]. 

Proof. According to Lemma f3.lt the function m_ (x, A) = e lXx f- (x, A) satisfies the Volterra equation 



(20) m-(x,\) = l+ K(x,y,X)m^(y,X)dy 

J — oo 

where K(x,y,X) := gjj (e 2 *^ - *) - l)V^ a {y). Note the bound 

\8%K(x,y,\)\ <C m \y-x\ m+1 \V^{y)\ 

and thus, Lemma fB.31 shows that the solution of Eq. (f20|) satisfies ||$™m_(-, A) || (— oo,a) < Cm for 
all A G [— Ao, Ao] since Ve t<T (y) decays exponentially as y — > — oo (Corollarv l2.1|) . The estimate for the 
derivative m'__(x, A) follows by differentiating Eq. (|2U|) and using the bounds for <9™m_(x, A). □ 

6.3. Zero energy resonances. Next we discuss the occurrence of resonances. First we give the 
precise definition of a zero energy resonance. 

Definition 6.1. We say that the operator Tie^ has a zero energy resonance if there exists a function 
/ G L°°(R) such that H i>a f = 0. 

Recall that the equation Hi^f = is equivalent to 

2M\ „ 2M , + 2Ma\ 
1 u + 5 1 ^— rw = 

where u(r) = f(x(r)) (see Eq. (jlip ). As already mentioned, solutions of this equation can be given 
in terms of special functions and therefore, we even know the behavior of f-(x, 0) for x — > 00. This 
is crucial to see whether the operator has a zero energy resonance or not. As the following 
lemma shows, no resonances occur for scalar perturbations. However, in the case of electromagnetic 
or gravitational perturbations one has to require £ > 1 (which we do anyway) or I > 2, respectively, 
in order to avoid resonances. 
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Lemma 6.3. Let (I, a) G N x {-3, 0, 1}\{(0, 0), (0, -3), (1, -3)}. Then the zero energy Jost 
solution has the asymptotic behavior f_(x,0) ~ cx f ~ +1 for x — > oo where c is a nonzero constant. 
In particular, there does not exist a function f G L°°(R) that satisfies Hi a f = 0- 

Proof. As already mentioned, the equation H£ a f = is equivalent to 

2M\ „ 2M , ( 1(1 + 1) 2Mo\ 

1 V 7TV + o 1 5— ) V = 



where v(r) := f{x(r)). Set u(z) := z"' 1+s )t)(2Mz) where s := \/l — °~ (note that s is the spin of 
the perturbing field, i.e., s G {0, 1, 2}). Then the above equation is equivalent to 

z(l — z)u + [c — (a + b + l)^]n — abu = 

where a = — £ + s, b = £ + 1 + s and c = 1 + 2s. This is the hypergeometric differential equation 
and we have the solution u(z) = 2 Fi(a, b; c; z) = %F\ {—I + s, £ + 1 + s; 1 + 2s; z) which reduces to 
a polynomial of order £ — s provided that £ — s G No (see p~]) and this is equivalent to (£, a) G" 
{(0,0), (0, —3), (1, —3)}. The Frobenius indices for the hypergeometric differential equation at the 
regular singular point z = 1 (which corresponds to r = 2M and hence, x — > —00) are (0, c — a — b) = 
(0,0) (see p]) which shows that u(l) 7^ and hence, there exists a nonzero constant cq such that 

f-(x, 0) = cor(x) 1+s u(^j-). Since u is a polynomial of order I — s, we obtain /_(x, 0) ~ cia/ +1 for 
x — > 00 by Lemma 12 . 1 1 where ci is a nonzero constant. □ 

Remark 6.1. Note that Lemma 16.31 is sharp in the sense that the operator Hi a does indeed have 
zero energy resonances if (£,cr) G {(0, 0), (0, — 3), (1, — 3)}. The resonance functions fi^ a are given 
by /o,o(z) = 1, /o-3(20 = 1 - 1(7) and /i,-3(sc) = as can be checked immediately. 

Remark 6.2. In the scalar case {a = 1), the hypergeometric function in the proof of Lemma 
reduces to the Legendre polynomial Pi and we have 

J V ' ; 2M V M 

6.4. Matching with /_(-,A). The above results are sufficient to match the Jost solution /„(-,A) 
to the solutions Uj(-,X) obtained in Lemma 14.31 by perturbing in energy. In what follows we will 
always assume that we are in the nonresonant regime, i.e., 

(£,a) £{(0,0), (0,-3), (1,-3)}. 
Lemma 6.4. The Wronskians cj (A) := W (/_(-, A), Uj(-, A)) for j = 0,1 are of the form^ 

cj(\) = 2i (l)+i0 2i+1 (\) 
for A G (0, Ao) where Xq > is a constant. Furthermore, we have cj~(0) ^ 0. 

Proof. According to Lemma l6.2( f—(x, A) is smooth in A around A = and by definition we have 
f-(x, A) = f-(x, —A) for A G 1. In particular, this implies Re/_(xo, A) = O^(l) and Im/_(xo, A) = 
02£+i{X) where xq > is the constant from Lemma 14.31 Repeating these arguments for the 
derivative f'_(x,X), we similarly obtain Ref'_(xo, X) = O2HI) and lmf'_(xo, X) = 02£+i(X). By 
construction (cf. Lemma [43]), we have uo(xo,X) = uo(xo), u' (xq,X) = u' (xq). Combining this 
with Corollary 14.21 we obtain Uj(xo, A) = O2KI) and u'j(xo, A) = (^(l) for j = 0, 1. This shows 

cj(X) = 2e (l)(0 2 i(l)+i0 2t+1 (X)) = 2 t(l) + i0 2i+1 (X). 

Suppose c^~(0) = W(/_(-, 0), ui(-)) = 0. This is equivalent to /_(-,0) and u\ being linearly depen- 
dent which implies that f-(x, 0) = 0(x~ i ) for x — > 00, a contradiction to Lemma 16.31 □ 



2 See Definition gT] 
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7. The spectral measure at zero energy 
Recall that we are interested in estimating the integral 

Acos(iA)Im [G e>a (x,x',X)] (x)- a (x')- a d\ 

and thus, we have to study the expressions 

/_(x',A)/+(x,A) 



Im- 



W(f.(;X),f + (;X)) 

In this section we obtain estimates for A — > 0+ and, as will be clear afterwards, the decay properties 
for solutions of the Regge- Wheeler equation are completely determined by this asymptotic behavior. 
In view of this, the following lemma is in fact the central result of our work. 

We have to consider different ranges of x, x' and A separately and we start with estimates for 
|A|, |xA| and |x'A| small which turns out to be the most important case. For all A £ (0, Ao) with a 
sufficiently small constant Ao > 0, we have the representation 

f±(x,X) = -cf(X)u (x,X) + Cq(A)ui(x, A) 

where c fW = W(f±,Uj)(X) (see Lemmas 16.11 and 16. 4p . Note carefully the slightly inconvenient 
fact that Cq (A) is the coefficient of ui(-, X) and not uo(-, A)! It follows that 

W(f., U)(X) = c7(A)c + (A) - Co -(A)c+(A). 

We abbreviate 

a m -t [ C 7( A )4(A) 
jki) [cr(A) C +(A)-c (A)c+(A)_ 
and, since Uj(-, A) are real-valued, we have to study expressions of the form 

A)o(A)«i(a;, A)tti(x', A), A 10 (X)u (x, Xju^x', A), etc. 
Lemma 7. 1 . The function Ajk is of the form 

Aqo(X) = 2m (A) and A jk (X) = 0(X 2i+1 ) if j + k > 1 
for X £ (0, Ao) where Ao > is a sufficiently small constant and the O-term behaves like a symbol. 

Proof. We have to distinguish four cases. 
(1) For Aqq we write 

%(A) 
Aw (A) = Im ^ 



l-f^(A) 
According to Lemma 16.41 we have 



£p_ fA) = c (A)c7(A) = Q 2£ (l) + iQ 2£+1 (X) 
q kr(A)| 2 |c7(A)P 

and, since |c7(A)| 2 = 2 e(l) + 2e+2 (X 2 ) = 2t (l) as well as |c7/(0)| / 0, we infer %(A) = 

c i 

2 e(l) +i0 2 £ + i(X) with the help of Lemma IA. 11 Furthermore, Lemma [6. II and Lemma |A. II 
yield 

£U A ) = _^ A 2^i (1 + 0(A£) + i0 {X- 2ie )) = 0(A 2£+1 - 2fe ) + iO(X 2e+1 ) 

Co a 



^From now on we write W(f±, Uj)(X) instead of W(f±(-, A), Uj(-, A)). 
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where the O-terms are of symbol type. Applying Lemma IA.1I again, we conclude 
A m {\) - 1m ________ 

= <WA) + 0(A 2 * +1 ) 

for a sufficiently small e and the O-term is of symbol type. 
(2) For A\q we use the representation 

Aio(A) = Im — -r = Im : 



c, ~cF^ ' 



_ l + Q(A^-^) + zO(A^ 1 ) _ 2m 
" Im i + 0(A 2£ + i-2 fe) " °( A ) 

where all O-terms are of symbol type (use Lemma IA.1|) . 
(3) In order to estimate Aq\ first note that 

W(f-,f+)(X) = cr(A)c+(A) - Co (A)c+(A) = ia (3 e c\- £ (l + 0(A £ ) + *0(A £ ^+ 2 ))) 

as follows straightforward from Lemmas 16.11 and 16.41 where c is a nonzero real constant and 
all O-terms are of symbol type. However, this implies 

„ m y ^^- c \ 2e+1 (0(l)+^0(\))(l + 0(\c) + iO(\-^)) 
Aoi(A) = Im— y 



1 + 0(A e ) + iO(A £ (^+ 2 ) 

Im 



i^X 2e+1 (0(l) + 0(A E ) + iO(A" 2fe ))(l + 0(A £ ) + iO(\< 2e+2 h) 



1 + 0(A £ 



Q(A 2m ) _ nfx z£+l\ 

1 + 0(A £ ) " ° (A } 



where all O-terms are of symbol type (see Lemma lA.ip . 
(4) Finally, for An we proceed exactly as above and obtain 

z^A^(Q(l) + ^(A))(l + 0(A £ ) + .0(A-^)) 
All(A) " Im l + 0(A £ ) + *0(A^+2)) "°( A ) 

where the O-term behaves like a symbol. 

□ 

Remark 7.1. The fact that j4qo(A) is somewhat exceptional is a direct consequence of the asymmetric 
decay properties of the Regge- Wheeler potential. This phenomenon is not present in |52j or |53| . 



8. Oscillatory integral estimates for small energies 

In this section we obtain bounds for the oscillatory integrals that describe the time evolution 
of solutions to the Regge-Wheeler equation. We distinguish different regimes, depending on the 
ranges of x, x' and A and in this section we only consider the case |A| small. As already mentioned, 
the most important contribution comes from the regime |xA| and \x'X\ small which yields the decay 
rates stated in Theorem 12.11 The remaining cases can be treated very similar to [53J, however, for 
the sake of completeness we give explicit proofs for all of the following statements. 
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8.1. Estimates for \xX\ and \x'X\ small. We will need the following result on oscillatory integrals. 

Lemma 8.1. For an N G No let w(A) = 02at+i(A) and suppose there exists a constant Ao > such 
that co(X) = for all A > Aq. Then we have the estimates 



f 

Jo 



Xcos(tX)ui(X)dX 

sin(iA)cj(A)dA 
for all t > where C{oo) can be estimated as 



< C(u)(t)-( 2N+ V 
<C( w )(t)- (2W+2) 



C(u) < Cmax |M|l°°(o,oo)> sup|A 



j^N+l+j) 



A>0 



(A)| :j = 1,2,3 



for an absolute constant C > 0. 



Proof. We only prove the sine estimate since the proof for the cosine estimate is completely analo- 
gous. It suffices to consider t > 1. (2N + 2)-fold integration by parts yields 

poo -i poo 

J sm(tX)u(X)dX = J sm(tX)J 2N+2 \X)dX 

since the boundary terms vanish thanks to oj^ 2m \0) = for m < N and the fact that oj(X) = for 
all A > Aq. Thus, it suffices to show that 



poo 

/ sm(tX)J 2N+2 \X)dX 
Jo 



< C 



for a constant C independent of t. Let x be a smooth cut-off satisfying x(A) = 1 for < A < \ 
and x(A) = for A > 2. Then we have 



poo poo 

\ sm(tX)J 2N+2) (X) X (tX)dX = / sm(r])J 
Jo Jo 



t 



10 JO 

for all t > 1. Furthermore, by an additional integration by parts we obtain 



sin(?7) 



dri<l 



f'OO 1 f'OQ 

/ sm(tX)J 2N+2 \X)[l - X {tX)]dX < - / cos(tX)J 2N+3 \X)[l- X (tX)}dX 
Jo t Jo 

+ - / cos(tX)J 2N+2 HX)t x '(tX)dX 
t Jo 

where the boundary term vanishes thanks to the cut-off and ti/ 2iV+2 )(A) = for all A > Ao- 
However, we have 



1 r 00 

- / cos(tAV( 2JV+3 )(A)[l-x(tA)]dA 
t Jo 



~ t 2 



t z 



cos(r]) — [1 -x(v)} 
rj z 



as well as 



1 f°° 

- / cos(iA)aj 
t Jo 



(27V+2), / 



tx'(tX)dX 



/ co S ( V )J 2N ^ (2) X '(V)^ 



< 



f 

Jo 



COS 7] 



since supp(x') C [5,2]. 



□ 
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Now we are ready to prove the first oscillatory integral estimate, valid for small A and |xA| < <5, 
|x'A| < 5 where 5 > is sufficiently small. In what follows we denote by xs a smooth cut-off 
function supported in a ^-neighborhood of the origin, i.e., 



Xs(x) 



1 if \x\ < | 



if \x\ > 6 ' 

Lemma 8.2. Let a > 21 + 1 and 5 > be sufficiently small. Then we have the estimates 



sup 



and 



sup 



A cos(tA)Im 



sin(iA)Im 



f.(x',X)f + (x,X) 

/-(x',A)/ + (x,A) 
W(/_(-,A),/+(-,A)) 



{x)- a (x>)- a X5WX5(xX)X5(x'\)d\ 



(x)- a (x')- a X8(\)X8(x\)x8(x'\)d\ 



<(t) 



-(2£+3) 



<(t) 



-(2£+2) 



for all t > 0. 



Proof. We set w(x, x', A) := Aqq(X)ui(x, X)ui(x' , A)(x) a (x') a xs(X)x5(xX)x5(x' X). According to 
Lemmas 17.11 and 14.51 we have w(x,x',A) = 02£+\(X) for fixed x,x'. Combining Lemmas 14.31 HIS 
andSSl we obtain \df +m ui(x, A)| < C m {x) 2£+l X~ m for m G N , A G (0,6) and x G [SX' 1 , SX^ 1 }. 
This implies 

\df+Mx,x',x)\<(x)^-«(xT +1 - a x-\ 

and, analogously, 

\2t+l-a /„l\2t+l-a \-2 



for all A > and x, x' G 



sup 



X cos(iA)u>(x, x , X)dX 



d 2 x e+3 [oj{x, x ',x)]\ < (x) 2t+ '- a (x'y zt+i - a x- 

\di^Hx, X ',x)]\ < (xr +i -°(x'r +i - a x- 3 

Finally, w(A) = for A > 5. Thus, Lemma 18 . 1 1 yields 
< (t)-(^+ 3 ) and sup 



x,x'& 



sin(iA)w(x, x' , A)dA 



(t) 



-(21+2) 



for all t > 0. For the remaining cases assume j ' + k > 1 and set 

u(x,x', A) := A jk (X)uj>(x, X)u k >(x', X){x)~ a {x')~ ol xs(^)X5(xX)xs(x'X). 

According to Lemmas 17.11 and 14.31 we have oj(x,x',X) = 0(X 2e+1 )0((xY +1 )0((x'Y +1 ) where the 
O-terms behave like symbols (use Proposition 14.11 and Lemma |4.6[) . In particular, oj(x, x',X) = 
C>2£+i(A) for fixed x,x' and 

\d 2 x i+2 co(x,x',X)\ < {xY +1 - a {x'Y +l - a X-\ 



as well as 



|df +3 Mx,*',A)]| < (xY +1 - a (x'Y +1 - Q x- 2 , 
\d 2 Y + *nx,x>,x)]\<(xY +1 - a (x'Y +1 - a x~ 3 



□ 



for all A > and x, x' G K. Thus, as before, applying Lemma 18.11 yields the claim 

Remark 8.1. Obviously, by performing fewer integrations by parts (cf. the proof of Lemma f8.ip . one 
may obtain weaker decay bounds (in t). By doing so, however, one can relax the decay requirements 
(in x) of the data, that is, the index a in Lemma 18.21 can be chosen smaller. To be more precise, 
one obtains the additional bounds, valid for a G N, a < 21 + 1, 

f-(x',X)f + (x,X) 



sup 



A cos(iA)Im 



W(/_(-,A),/+(-,A)) 



(x)- a (x')- a X5Wx 5 (xX)x S (x'X)dX 



(t) 



-a-2 



2U 



and 



sup 



for all t > 0. 



sin(tA)Im 



f-(x>,\)f + (x,\) 

W(f_(;X)J + (;X)) 



(x)- a (x>y a X8WX5(x\)x5(x'\)d\ 



-a-1 



Remark 8.2. Note that the sine estimate from Lemma 18.21 is the main obstacle to proving better 
decay. The remaining oscillatory estimates of Sec. [8] and Sec. [9] below are consistent with faster 
decay. 

8.2. Estimates for \xX\ and \x'X\ large. For the remaining small energy contributions it is useful 
to note that, for A G R, f±(x,—X) = f±(x,X) by definition of the Jost solutions. This implies 
Ge !(T (x, x' , — A) = Gi :(7 (x, x', A) and hence, the real part Re [G£ )CT (x, x', A)] is an even function of A 
whereas the imaginary part Im [Ge t0 -(x, x' , A)] is odd. Thus, we have 

/ A cos(iA)Im [Ge :a (x, x', A)] dX = - / Xcos(tX)Gi iCT (x,x' , A)<iA 
Jo "2 J R 

and similarly for the sine evolution. This shows that we can replace the imaginary part of 
Gp ja (x, x' , A) by G^ jtT (x, x' , A) itself in the oscillatory integrals and change the domain of integration 
from A > to A 6 1. Furthermore, recall 



(21) 



W(/_,/ + )(A) = cr(A)c + (A) -Cq (A)c+(A) = cA^(l + O c (A e )) 



for a nonzero constant c where the O-term is of symbol type. This has been shown in the proof of 
Lemma 17.11 

In order to deal with terms that involve f-(x', A) for x' > and f+(x, A) for x < we have to 
consider reflection and transmission coefficients. For A 7^ 0, the functions /+(-,A) and /+(-,A) are 
linearly independent which shows that there exist coefficients a(A) and 6(A) □ such that f-(x, A) = 
a(X)f + (x, A)+6(A)/ + (x, A). This representation implies |6(A)| 2 — |a(A)| 2 = 1 (cf. the proof of Lemma 
K2b and thus, /+(x,A) = — a(A)/_ (x, A) + 6(A) /_(», A). Furthermore, we have W(/_,/+)(A) = 



1 

2iA' 



b(X)W(f + , /+)(A) = 2iA6(A) which is equivalent to 

6(A) 
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W(/-,/+)(A) 

Similarly, we obtain W(/_,/ + )(A) = — 2iAa(A) and therefore, 

o(A) W(/_,7^(A) 



W(/_,/+)(A) 2iAW(/_,/+)(A)' 



However, from Lemma HTT1 and W(/_, /+)(A) = c : (A)cq (A) — c (A)q (A) it follows that 

^(/-,/ + )(A) =C + 0c(A) 
where c is a nonzero constant and the O-term is of symbol type and hence, 

a(X) 



(23) 



w(f-,f+)(X) 



cA- 1 (l + O c (A £ )). 



In order to avoid confusion, we remark that the coefficients a and b axe not exactly the same as A and B in the 
proof of Lemma 13.21 but they are related by a(A) = A(\ 2 ) and 6(A) = B(X 2 ). 
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Lemma 8.3. Let a E No and 5 > sufficiently small. Then we have the estimates 



sup 

x,x'eH 



Acos(tA) ^^4^ (x)- a (x')- a X s(\)(l - Xs(xX))(l - X s(x'X))dX 



and 



sup 



sin(tA) 



W(/_,/+)(A) 

/_(x',A)/ + (x,A) 
W(f.,f + )(X) 



< /+\-a 



(x)- a (x')- a Xs(X)(l - Xs(xX))(l - xs(x'X))dX 



for all t > 0. 



Proof. Let |A| < 5, x > 0, xf < and |Ax|, \Xx'\ > f. We set m±{x, A) := e TiAa 7 ± (x, A). According 
to Lemma 16.21 we have the bound |<9™m_(x,A)| < C m for all x < and m G No- Furthermore, 
since A is small, we have 



m+(x, A) = e" 



iXx 



x, A) = (1 + O c ((Xx)- L ))(l + h(Xx, A)) 



by Lemma 15.11 and the asymptotics of the Hankel function where the O-term is of symbol type. 
By Lemma l5?2l and the chain rule, we have the estimate \d x 1 b£(Xx, A)| < C m (x} m since |A| _1 < (x) 
and this implies \d™m+(x, X)\ < C m (x) m for all m 6 No- Set 

w(x,x',A) := A ^4yW(^,A) ^ (A)(1 _ w(a?A))(1 _ w(a;/A)) 

Then we have 

|flS» w (x,a/,A)|<<x) m <a'r 
for all m G No since W(/_, /+)(A) is of symbol type (cf. Eq. (J2Tj) ). We have to estimate the integral 

3 iA(±t+x-x') w ( a . j x / ; A)(x)- a (x')- a dA. 



If | ± t + x — x'\ > \t we integrate by parts a-times to obtain 



< 



and, if | ± f + x — x'| < |f, we have (x) °(x') Q < (t) a as f — )• oo and thus, 



iA(±t+a;-a/) 



w(x,ar / ,A)(x>- a (a/>- a dA 



<*>" 



If x' > or x < we use the representations m_(x',A) = a(X)e 2tXx m+(x', A) + 6(A)m + (x',A) 
or m+(x,A) = —a(X)e~ 2iXx m_ (x, A) + 6(A)m_(x, A) and with the help of Eqs. (j22j) and (J23J) the 
corresponding integrals can be estimated as above. 

For the sine evolution note that we are missing one A and thus, for instance, if 

, a(A)m + (x',A)m + (x,A) , 
w(x, x , A) := — - Xi(A)(l - x*(xA))(l - Xa(z A)) 

we have |w(x,x',A)| < | A | ^ 1 (cf. Eq. ([21])) which, in the case t = 0, only yields the weaker bound 
|3^w(x,x',A)| < (x) m+1 (x') m . □ 
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8.3. Estimates for |xA| small and \x'X\ large. The next regime to be considered is \x\\ < 5 and 
A| > f . We distinguish the two cases x' < and x' > 0. 

Lemma 8.4. Let a G N and 5 > sufficiently small. Then we have the estimates 

J-(x',\)U(x,\) 



and 



sup 

x&S.,x'<0 



sup 

xeM,a:'<0 



A cos(iA) 



/ sin(tA) 
Jr 



w(f.,f + )(\) 

f-(x>,\)f + (x,\) 

w(f.j + )(x) 



■(x}- a (x'y a X s(X)X6(xX)(l - Xs(x'X))dX 



< (t)- a 



(x)- a (x')- a xs(X)x6(xX)(l - X s(x'X))dX 



< (t)- a+1 



for all t > 0. 



Proof. Let x G R, x' < 0, |A| < 5, \Xx\ < 5 and |Ax'| > 2- As before, we write f±(x,X) = 
e TtXx m±(x, A) and recall the bound |(9™m_(x', A)| < C m for all m G No (Lemma 16. 2|) . For m + we 
use the representation 

m+(x, A) = e~ iAa; (-c ] f '(A)no(x, A) + c^{X)ui(x, A)) 



where we extend Cj (A) and Uj(x,X) to negative A according to cj~ (— A) = c+ (A) and Uj(x,— A) 
Uj(x,X). Applying Proposition 14.11 and Lemmas 14.61 16.ll we obtain the bounds 



\d^ct(X)u (x,X)\ < (x) e+l \Xf- m < (x)\X\- m < (xT\M~ l 



and 



m < 



\d?4(x) ui ( x ,x)\ < (x)\x\- e - m < (x'riM-*- 1 

for m G No which implies \d™m + (x, A)| < {x') m \X\~ l ~ l . We also have Id^xsixX)] < (x) 
\X\~ rn < (x') m and, putting all this together, we arrive at 

\d?u(x,x',X)\<(xT 

where 

t i \\ ,m-{x',X)m + (x,X) , 
u{x, x,X):= X w X5WX5{xX){l - xs(x A)), 

see also Eq. f)21 1) . Thus, the claim follows by appropriate integration by parts as in the proof of 
Lemma 18.31 The proof for the sine evolution goes along the same lines but one loses one power of 
A. 

Lemma 8.5. Let a £ N and 5 > sufficiently small. Then we have the estimates 

f (x',A)/ + (x,A) 



□ 



and 



sup 

x&R,x'>0 



sup 

x£R,x'>0 



/ Acos(tA) 
Jr 



sin(tA) 



w(f.,f + )(X) 

f_(x',X)f + (x,X) 
W(f.J + )(X) 



■(x)- a (x')- a xs(X)xs(xX)(l - X s(x'X))dX 



<(*)' 



(x)- a (x')- a X8(X)x8(xX)(l - xs(xX))dX 



(t) 



-a+l 



for all t > 0. 



Proof. Let x G E, x' > 0, |A| < 5, \xX\ < 5 and \x'X\ > |. As always, we write f±(x,X) 
e^ lXx m±(x, A). Again, by Lemmas 14.31 16. l l 14.61 and the representation 



we obtain the bound 



m+(x, A) = e- lAx (-cj(X)u (x, A) + c£(X) Ul (x, A)) 
\m + (x,X)\ <(x) e+1 \Xf + (x)\X\-*<(x)\X\-t 
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and, by using the symbol behavior of the involved terms (see Proposition 14. II and Lemma [6. ip . this 
implies 

\d?m+{x,\)\ < (x)|A|-^ m < (x)(xT\M- e - 
For m_(x',A) we use reflection and transmission coefficients, i.e., 

m-(x', A) = a{X)e 2iXx ' m + {x' , A) + b(X)m + (x', A) 

and, from the proof of Lemma [8.3l we have |(9™m+(x', A)| < {x') m . Hence, Eqs. (j23|) and (|22[) show 
that 

\d?m-(x>, a)| < iai- 1 -™ + (xtim- 1 < (x'riM- 1 - 

Setting 

u){x, x , A) := A w y / )(A) ^(A)x<5(^A)(l - A)) 

the above estimates and Eq. (|2ip imply |<9™w(x, x', A)| < (x){x') m and the cosine estimate follows 
by appropriate integration by parts as in the proof of Lemma 18.31 For the sine estimate we set 

u(x, x , A) := ^— XsWxs(xX){l - xs(x A)) 

and we only have the weaker bound \d™u(x, x', A)| < (x)(x') m |A| _1 < {x){x') m+1 . □ 

8.4. Estimates for |x'A| small and \x\\ large. Due to the asymmetric decay of the Regge- 
Wheeler potential, this case is slightly different from the above considered \x\\ small and |x'A| 
large. Thus, it has to be studied separately and does not follow from symmetry arguments as in 
J. Similar as above, we distinguish x > and x < 0. 



Lemma 8.6. Let a £ N and 5 > sufficiently small. Then we have the estimates 
Acos(tA) / -5' A) {^' A) (x)- a (x / )-^(A)(l - X s(xX))xs(x'X)d\ 



and 



sup 

x>0,x'& 



sup 

x>0,x'eR 



sin(iA) 



f-{x',X)f+(x,X) 
W(f.J + )(X) 



<*>■ 



(x)- a (xy a xs(X)(l - X s(xX))xs(x'X)dX 



< (t)- a+1 



for all t > 0. 



Proof. Let x > 0, x' 6 M, |A| < <5, |x'A| < <5 and |xA| > |. As in the proof of Lemma 18.31 we have 
the bounds \d™m+(x, A)| < (x) m . For m_(x', A) we use the representation 

m_(x', A) = e iA:r '(-cr (A)uo(x', A) + c (A)ui(x', A)) 
since |x'A| is small. Lemmas 14.31 14.61 and 16.41 imply the bound 

\m-&,\)\<W* 1 + (J)<\\\- t - 1 

and, by using the symbol behavior (see Proposition 14.11 and Lemmas 14.41 16. 4p and (x') < |A| _1 , we 
infer 

\d?m-(x',\)\ < lAI^- 1 — < (x)-lAI^- 1 . 

We set 

( i w \ m -( x ', X)m + (x, A) , 
u)(x, x , A) := A — — — — — — - — x<j(A)(1 - Xs{xX))x5{x A) 

and the above bounds as well as Eq. (|2Tj) imply the estimate \d^u:(x,x',X)\ < (x) m . Thus, as 
before, the claim follows by appropriate integration by parts. □ 
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Lemma 8.7. Let a G N and 5 > sufficiently small. Then we have the estimates 



and 



sup 

x<0,x'& 



sup 

x<0,x'£R 



\co S (t\) f -(f/ f X) { + [^\ x}- a (x'r a xs(\)(l - X s(x\))xs(x'\)d\ 



sin(iA) 



W(/-,/+)(A) 
f-(x',\)U(x,\), . a . , 



W(f.J + )(X) 



(x)~ a (x')- a Xs (X)(l - X s(xX))x5(x'X)dX 



< (t)- a+1 



for all t > 0. 



Proof. Let x < 0, x' G R, | A| < (5, |x'A| < 5 and |xA| > I. Like in the proof of Lemma IHTH we have 



\d^m.(x',X)\<{x') e+1 \X\- m <(x')\X\- 



-l-m 



for all m G No- Since x < we use reflection and transmission coefficients to obtain the represen- 
tation 

m+ (x, A) = -a(X)e~ 2iXx m^(x, A) + 6(A)m_(x, A) 



which immediately implies the bound \m + (x,A)| < (Aj- 1 by Eqs. d22D and LemmaE3 Thus, 
from the symbol behavior of a(A), 6(A) and \d™ni-(x, A)| < C m (Lemma 16. 2p . we infer 

\dTm + (x,X)\ ZixriM" 1 
for all m G No since |A| _1 < (x). Thus, Eq. (|2Tj) implies 

\8?u(x,a/,X)\ <(x) m (x') 

for all m G No where 

, u m-(x',X)m + (x,X) , u . , , » 

w(ar, x , A) := A ^. Xa(A)(l - x<5(zA))x5(x A). 

As a consequence, by appropriate integration by parts (cf. the proof of Lemma f8T3|) . we obtain 



JX(±t+x-x') 



uj(x,x',X){x)- a {x')- a dX 



(ty 



provided that a > 1. For the sine evolution we set 



and use the bounds 



to obtain 

i^(^^,A)i<(xr +i (x') 

for all m G Nq. Hence, as before, the claim follows from 



J\(±t+x-x') 



u{x,x',X)(x)- a {x')- a dX 



which can be obtained by appropriate integration by parts similar to the proof of Lemma 18.31 □ 
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9. Oscillatory integral estimates for large energies 

9.1. The Jost solutions at large energies. In order to estimate the contributions from large 
energies, we need the behavior of the Jost solutions for A — > oo. As usual, we write m±(x,X) = 
e^ iXx f±(x,\). 

Lemma 9.1. Let Ao > 0. Then, for k,m G No, the function m + (-, A) satisfies the estimates 

\d*d?(m+(x,\) ~ 1)1 < C k , m {x)- l - k \- l - m 
for all A > Ao and all x > 0. The same bounds hold for m_(x, A) if x < 0. 

Proof. As already discussed (see Lemma l3.ip . the function m+(-,A) satisfies the Volterra equation 
m+(x, A) = 1 + — J (e 2iA ^) - 1) V^(y)m+(y, X)dy 

poo 

= 1+ (e 2 "> - 1) 
Jo 

and thus, the Lemma is obviously true for k = m = 0. Let n : No X No — > No denote the 
bijection from Proposition 14.11 Fix (k, m) 6 No x No and suppose the assertion is true for all 
(J, I) with n(j,l) < n(k,m). We need to show that this implies the claim for (k',m') where 
n(k',m') = n(k,m) + 1. There are two possibilities: Either (k',m') = (m + 1,0) (if k = 0) or 
(£/, m!) = {k — l,m + 1). In the former case we have 

dr l m + (x, A) = -L ( m j ) I " 1) dr 1 -^ (I + x) dlm + Q + x, X) dr, 

+ 2^ i°° ^ " x) V ^ (x + x ) (x + x ' A ) dr? 

= Ocdxy'-^X- 1 ) + (e 2lA ^ - l) V e>rT (y)d™ +1 m + (y, X)dy 

by assumption and thus, the estimate from Lemma lF3.ll yields |<9™ +1 m + (x, A)| < (x) _1_ ( m+1 )A _1 . 
For the latter case observe that 

and, more generally 0, 

d[dim + (1 +X ,X)<(1 + ,) ' ' ' ± iX ~l + • • • + (2 + A— < (| + x)" W A- 



r? \ m+ (% + x > A ) 
A / 2iA 2 



drj 



^One may apply Faa di Bruno's formula (cf. Lemma [A. 1|) to obtain a completely explicit expression for the higher 
A-derivatives. 
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for all (j, I) with n(j, I) < n(k, m) and (j, I) ^ (0, 0) by assumption. This shows that 



E 

j=0 



m + 1 

j 



2i V 



+ / (e 2 "> - 1) d k x 



k-l 



1) dt 1 



2iA 2 



5 



m+l 



2iA 2 



-+x,X 



rriA 



+ x, A 



r/r/ 



(/?/ 



Oc((ar)- 1 - (fc - 1) A- 1 -( m+1 )) 



+ 



2iA 



and Lemma |B , 1 1 yields the claim. The proof for m_ is (mutatis mutandis) identical. 
Corollary 9.1. Let Xq > 0. Then the Wronskian W(/_,/ + )(A) has the behavior 

-(l + O c (A- 1 )) 



□ 



W(f-,f+)(X) 2iX' 
for all A > Aq where the O-term behaves like a symbol. 

Proof. With m±(x,X) = f±(x, A) we have 

W(f-, /+)(A) = 2iAm_(0, A)m+(0, A) + Oc^ 1 ) = 2iX(l + Oc(A" 1 )) 

by Lemma 19 . 1 1 where the O-term is of symbol type. Thus, the claim follows from Lemma lA.ll □ 

Before proceeding to the final oscillatory integral estimate, we need the large A behavior of the 
reflection and transmission coefficients, i.e., the coefficients a(A) and 6(A) satisfying f-(x, A) = 
a(X)f+(x, A) + b(X)f+(x, A). The behavior of b is given by Eq. (|22j) . By Lemma 19. II we have 

W(f-,U)(\) =m_(0,A)m' + (0,A) - m'_(0, A)m+(0, A) = Oc^ 1 ) 

where the O-term behaves like a symbol and therefore, 

a(X) W(f„,U)(X) 



(24) 



O c (A" 



W(f-,f + )(X) 2iXW(f.J + )(X) 

We also remark that, by symmetry, the above considerations extend to large negative A. Now we 
are ready to prove the final oscillatory integral estimate. 

9.2. The cosine estimate. We distinguish between the cosine and the sine estimate since in the 
former case we obtain a bound involving the derivative of the data. 

Proposition 9.1. Let a E No and 5 > sufficiently small. Then we have the estimate 

i-N 



sup 



lim 

Af-5-oo 



R J-N 



Xe ±ux G^ a {x, x', A)(l - X s{X))(x)- a (x')- a 4){x')dx'dX 



< 



V)I + I<K*')I)^' 



for all t > and any 4> £ <S(M). 

Proof. We split the integral according to 

/ / ... dXdx' 
Jr J-n 



N 



, dXdx + 



-oo J-N 



oo pN 
N 



. dXdx 
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and only consider the first summand since the proof for the second one is completely analogous. In 
the domain x' < x, which we study now, the Green's function is given by 



Ge a (x,x',X) 



f.(x',X)f + (x,X) 



We distinguish between x < and x > and start with x < 0. Using reflection and transmission 
coefficients we obtain 



/_(z',A)/ + (x,A) = e- iXx 'm-(x',X) —a(X)e~ iXx rri-(x, A) + b(X)e iXx m-(x, X) 

= -e~ iX{ - x+x ' ) 'aJX)m- (x', A)m_ (x, A) + e iX ^ x - x %{X)m^ (x' , A)m_ (x, A) , 

and consider each term separately. We define 

, Xb(X)(l- Xs (X)) , — r — 
ui{x,x,X) := w y f )(\) rn-{x ,X)m^{x,X) 

and by Lemma 19. II and Corollary 19.11 as well as Eqs. (|22p . (|24|) . we obtain the estimates 

\d?u(x,x',\)\ <c m \xr m 

for all |A| > | and x,x' < (recall that 1 — Xs(X) = 1 for |A| > 5). Note that, by Fubini, we can 
freely interchange the order of integration and thus, integration by parts with respect to x' yields 

(25) 



x pN 



z iX{±t+x - x,) Lu{x,x',X){x)- a {x'}- a (j)(x')dXdx' = -<f>(x){x) 



-2a 



oo J-N 



N ±i\t 



N 



iX 



lu(x, x, X)dX 



+ 



x pN i\(±t+x-x') 



oo J—N 



iX 



-{x)- a d x , [uj{x,x',X){x')- a 0(x')] dXdx 1 . 



We first claim that 
(26) 

Indeed, set 



sup 

x<0 



e ±iXt X- 1 u(x,x,X)dX < 1. 



6(A)(l-x«(A)), 



u>{x,X) := A 1 u(x,x,X) = ^^y— ^^^m_(x, A)m_(x, A) 

and observe that u(x, X) = 2^(l + Oc(|A| -1 )) f° r l-^l — ^ oo by Eq. ([22]) and Lemma [97TT This shows 
that co(x, —A) = — Co[x, A) + Oc(|A|~ 2 ). Thus, we have 



f e ±iXt u(x, 



X)dX 



< 



< 



e ±iXt iv(x,X)+e TtM Cj(x,-X) 



=F*At; 



dX 



sm(Xt)Cj(x, X)dX 



+ 



e ±iXt O c (\X\-' z )dX 



sin(At)(A- 1 +Oc(|A|- 2 ))(iA 



+ 1<1 
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for all x < and this proves Eq. (|26p . Therefore, for N — > oo, we can estimate the first term in 
Eq. d25]) as 



4>(x)(x)- 2a / e ±lXt \- 1 oj{x,x,X)dX 



< sup\(p{x) 

x<0 



1 



)|- / e ±^[A-V^,A)] 



dX 



< 



(t)- a / {\<j>i(x')\ + \cf>(x')\)dx' 



for all t > 1 and x < by a-fold integration by parts and Sobolev embedding. By Eq. (|26p this 
inequality is in fact valid for all t > 0. 

For the second term we similarly claim that 



(27) 



sup 

2<0 J — OO 



fix' 



< 



(x)- a (xr a {W(x>)\ + \^x')\)dx'. 



Indeed, we have 



f jK±t+*-*')\-i( x )-*d x ,u(x,x',X){xY a <f)(x')dX <(x)- a (x')- a [ \X- x d x ,u{x,x',X)<l){x')\dX 
Jr Jr 

< (x)- a (x')- a W)\ 
for all x < since | X~ 1 d x tuj{x, x' , X)\ < |A|~ 2 by Lemma 19.11 Moreover, 



z iXi - ±t+x - x '^ X- 1 oj{x,x' ,X){x)~ a d x , [{x')- a (t>{x')} dX 



{x)-"(x')- a {W{x')\ + \<t>{x')\) 



for all x < which can be shown by exploiting exactly the same cancellation that led to Eq. (J26J). 
This proves Eq. (|27p . Note in particular that Eq. (|27p implies 



lim 

jV->oo 



x pN e i\(±t+x-x') 



— x -/-A ^A 

g iA(±t+x-x') 



/A 



(x>- Q cV [w(x,x',A)(x , )" Q ^(^)] rfAda/ 
(x>- Q cV [u(x,x',X)(x'y a 4>(x')] dXdx' 



by dominated convergence. Now we distinguish two cases. If | ± t + x — x'\ > |i, we integrate by 
parts a-times to obtain 

e i\(±t+x-x>) x -l( x ya dx/ [ u) ^ XjX ' j X){ X ')- a (l)(x')] dXdx' 

< | ± / + <■ - r'r" I I \{x)- a d*d% [X- 1 cj(x,x',X)(x')- a <t>(x , )]\dXdx' 



< 



(t)- a / \X\-( a+1 \x)- a \d xl [{x')- a Hx')}\dXdx' 
JrJ\\\>^ 

(t)~ a [ [W{x>)\ + W)\)dx> 



: : !S 



for all t > and all x < 0. If | ± t + x - x'\ < \t, we have < (t)~ a as t ->• oo and 

Eq. (f2T|) implies 



sup 

x<0 J — oo 



s iA(±t+x-x') A -l^-a 5x< [^^A)^')^^')] rfA 



dx 



< 



{t)~ a / [W{x')\ + \cj>{x')\)dx'. 



for all i > 0. The remaining cases are treated in a completely analogous fashion. Note that the 
terms involving the coefficient a(A) are even simpler due to the stronger decay given by Eq. (|24p . For 
terms that contain no reflection and transmission coefficients, use Corollary 19. II for the cancellation 
argument. □ 

9.3. The sine estimate. The sine estimate is slightly stronger since it does not require derivatives 
of the data. 

Corollary 9.2. Let a G No and 5 > sufficiently small. Then we have the estimate 



sup 



lim 

jV-s-oo 



N 



e ±ttx G e , a (x, x' , A) (1 - xs (A) ) (x) ~ a (x'^^dx'dX 

< (t)- a [ \cf>(x f )\dx' 



for all t > and any <fi G <S(R). 

Proof. Just repeat the arguments from the proof of Proposition 19.11 However, note that we are 
lacking one factor of A compared to Proposition 19.11 which makes the integration by parts with 
respect to x' unnecessary. This explains why no term containing <j)' appears on the right-hand side 
of the estimate. □ 



Appendix A. Symbol behavior 

Lemma A.l. Let Jcl and suppose f is smooth on L and satisfies |/(x)| < C < 1 for all x £ L. 
Then, for all x £ L, we have the estimate 



l + f 



(k) 



(x) 



< c *En|/ (i) ( a 

3=1 



for all k G N where the sum runs over all possible k-tuples (mi, ?7t2, • • • , rnk) G Nq satisfying 



Proof. This follows from the identity 



l + f 



(AO 



E 



"mi ,rrt2,...,rrifc 



n f/ w ) 



1 + / 



which is known as Faa di Bruno's formula (see e.g. [50], the explicit form of the coefficients 



is irrelevant for our purposes) and the fact that |(1 + /) 1 \ < 1 on /. 
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□ 



Appendix B. Volterra integral equations 

In this section we establish some well-known facts about Volterra integral equations which are 
frequently used throughout this work. 

Lemma B.l. Let a G R, g G L°°(a,oo) and suppose the integral kernel K satisfies 

/•CO 

\i := / sup \K(x, y)\dy < oo. 

J a x£(a,y) 

Then the Volterra equation 

/•oo 

f(x)=g(x) + K(x,y)f(y)dy 

J X 

has a unique solution f satisfying 

||/||l,°°(a,oo) < eM ||5lU°°(a,oo)- 

Proof. See e.g. |25] or [52]. □ 

The next lemma states differentiability properties of solutions of Volterra integral equations. 

Lemma B.2. If, in addition to the assumptions of Lemma \B.ll g G C°°(a, oo) and the kernel K 
is smooth in both variables on (a, oo) and satisfies 

oo 

\d k K(x,y)\dy < oo 

for any x > a and all k £ N then the solution f is smooth on (a,oo). Furthermore, the derivatives 
can be calculated by formal differentiation, i.e., 

1 poo 

f( k \x) = g( k \x) - ^(^/) (fc - 1 - j) (x) + / d k x K(x,y)f(y)dy 

3=0 Jx 



where Kj(x) := dlK(x,y)\ 



y=x- 



Proof. The claim follows from a straightforward application of Lebesgue's theorem on dominated 
convergence and an induction. □ 

The next lemma shows how the dependence of the kernel K on a parameter A carries over to the 
solution of the corresponding Volterra equation. 

Lemma B.3. Let I Cl be open and suppose 

/•oo 

/ sup \d^K(x,y,X)\dy < oo 

J a x£(a,y) 

as well as d™g(-, A) E L°°{a, oo) for all m G No and A G I. Then the Volterra equation 

/•CO 

f(x, A) = g(x, A) + / K(x, y, X)f(y, X)dy 



has a unique solution f(x, A) for all x > a and X £ I which is smooth in X. Furthermore, we have 
d™f(-, A) G L°°(a, oo) for all m G No and the derivatives are given by 

d?f(x,X) = d?g(x,X) + J2{™ ) J cV x K(x,y,X)d^f(y,X)dy. 



j=0 
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Proof. According to Lemma[Bj!l the solution / exists, is unique and satisfies ||/(-, A)||£ao( O 0O ) < oo. 
Now consider the integral equation 

poo 

(28) h(x, A) = d\g(x, A) + h(x, A) + / K(x,y, X)h(y, X)dy 

J X 

where 

/•oo 

h(x,X):= d\K(x,y,X)f(y,X)dy. 



We have \\h(-,X)\\ 

L°°(a, oo) < 00 an d thus, by Lemma [BJ4 Eq. (128j) has a unique solution h(-,X) £ 
L°°(a, 00) for all XE I. However, by dominated convergence we conclude 

f(x,X + v)-f(x,X) 



lim 

i/-s>0 



A) 







and hence, d\f exists and equals h. Existence of the higher derivatives follows by the Leibniz rule 
and an induction. □ 

We finally remark that all of the above Lemmas have counterparts for Volterra equations of the 
form 

f(x, A) = g(x, A) + f K(x, y, X)f(y, X)dy 

J a 

with almost identical proofs. 
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